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HYDRODYNAMICS OF LIQUID HELIUM II* 


C. C. Lin 


Department of Mathematics, 
Massachusetts Institute of Technology, 
Cambridge, Massachusetts 
(Received February 24, 1959) 


This note gives a brief account of a macro- 
scopic hydrodynamic theory of liquid helium I 
which can describe most of the experimental ob- 
servations on the thermal and mechanical be- 
havior of liquid helium below the A-point. The 
theory evolves from a one-fluid concept of liquid 
helium in which a long-range order exists in the 
momentum space. The transfer of long-range 
order (disorder) is then the superconductivity of 
heat. By assuming that one is dealing with a 
classical thermomechanical medium, it can be 
mathematically proved that there is a momentum 
transfer associated with the reversible energy 
transfer even when the momentum of the fluid is 
zero. This is a typical feature of the usual two- 
fluid theory; indeed, a two-fluid model can be 
constructed, in a purely mathematical manner, 
out of the present one-fluid theory. The basic 
equations, one governing the rate of change of 
momentum and the other governing the rate of 
change of the reversible heat flux vector, can 
then be reduced to the ordinary equations of 
Landau. There is, however, no longer any com- 
pelling reason to believe that the motion of the 
superfluid component must be irrotational. (See 
Lee and Yang.’) 

Accordingly, we propose to abandon the often- 
used basic assumption that the superfluid com- 
ponent of the flow is irrotational. (See Lifsic 





and Halatnikov.”) We next include certain dissi- 
pative terms in our equations of motion depend- 
ent on the gradients of physical quantities, in 
agreement with Lifsic and Halatnikov. Since the 
equations of motion now contain the second-order 
space derivatives of the superfluid component, 
we need a boundary condition for the tangential 
component of the superfluid velocity. Bearing in 
mind the requirement of covariance with respect 
to coordinate transformations, we propose the 
following general form for the boundary condi- 
tion: 


ov 7/®N =f@)4, (1) 


where q is the tangential component of the differ- 
ence between superfluid velocity and wall velocity, 
and 8V, 7/@N is the tangential component of the 
normal derivative of the superfluid velocity. 
Before a molecular theory can be developed, the 
form of f(@) must be determined by experiments. 
For small velocities, we assume 


f@)=a+6G+... . (2) 


From Andronikashvili’s disk pile experiment, 
we know that a must be negligibly small. Thus, 
the simplest form for f(q@) is 


f @) =Be. (3) 


With this simple assumption, we can derive the 
following results for experiments involving sim- 
ple geometry. 

(A) For rotating bucket experiments, the 
theory predicts uniform rotation of the super- 
fluid component, but at a speed lower than that 
of the normal component, because of the bound- 
ary condition (1). The ratio of angular momen- 
tum L actually acquired by helium II to the ratio 
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FIG. 1. Angular momentum of liquid helium II in a 


rotating bucket. 


L, for ordinary fluids is given by 
L/L, =s+x (1-s), 
w/w, = 2 [s/(1-s)?, (4) 
where x is the fraction of normal fluid present, 


w is the angular velocity of the bucket, and wy, 
is related to the molecular parameter 8 by 


wr =4Ba->. (5) 


In this formula, a is the radius of the bucket. 

(B) For isothermal flow through a tube (or 
channel), the amount of flow depends on the pres- 
sure gradient [ in the following manner: 


Mean velocity 7 =AdT + Bd“* Tr", (6) 


where A and B are molecular parameters, and d 
is the diameter of the tube (or width of the chan- 
nel). 

In Fig. 1, one set of experimental results of 
Walmsley and Lane*® is compared with Eq. (4). 
There is only one other set of data at the same 
temperature which is, however, not sufficiently 
accurate to test the validity of (5). In Fig. 2, the 
experimental results of Atkins* are compared 
with (6). One may notice that the theoretical and 
experimental results agree for tubes of various 
sizes. 

Preliminary analyses also show agreement of 
the theory with experimental results for fountain 
effect and for damping experiments at finite am- 
plitudes.° These will be discussed in a more 
complete paper. It also follows from the present 
concepts that second sound should propagate 
through rotating liquid helium II in essentially 
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FIG. 2. Velocity of flow of liquid helium II through 
tubes of various diameters. (I) d=4.40x 107? cm (Vv 


for short tube). (II) d=2.03x10~* cm (A for short tube), 


(II) d= 0. 815x107? cm. 


the same manner as in the nonrotating case—a 
fact observed by Hall and Vinen.® The detailed 
examination of attenuation of second sound in 
rotating helium is left for further investigation; 
preliminary considerations also show general 
agreement with their results. 

I should like to express my thanks to Professor 
C. N. Yang for many discussions during the de- 
velopment of the one-fluid concept of liquid he- 
lium II while both of us were at Madison during 
the summer of 1958. I wish also to thank Pro- 
fessor C. T. Lane for private communications 
and comments on the data shown in Fig. 1, and 
to thank Mr. Alfred Clark, Jr., for his help in 
the numerical analysis of the experimental data 
shown in Fig. 2. Thanks are also due to my 
colleagues L. Tisza, Kerson Huang, and Charles 
E. Chase for their interest and discussions. 
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MAGNETIC SUSCEPTIBILITY OF PHOTO- 
GENERATED CURRENT CARRIERS 
IN GERMANIUM 


J. O. Kessler and A. R. Moore 
RCA Laboratories, 
Princeton, New Jersey 
(Received February 4,. 1959) 


The experiments described in this note relate 
to measurements of the magnetic susceptibility, 
y, of free carriers in germanium by means of a 
method fundamentally different from those here- 
tofore employed. We observe a carrier y sev- 
eral hundred times bigger than that expected in 
the light of previous work reported in the recent 
literature. In the past, the carrier y has been 
derived by subtracting the y measured for a 
pure specimen from the x of a heavily doped one. 
In another approach, the temperature of an in- 
trinsic specimen has been raised so as to intro- 
duce many pairs in equilibrium. The carrier x 
was then obtained by subtracting the low-tem- 
perature, or “carrier-free” y from the total 
high-temperature x.” 

The results of these experiments fit fairly 
well the expectations derived from the band 
structure. Nevertheless, the validity of the sub- 
tractions is open to some general doubts.° 

The experiments reported here are based on a 
new experimental technique which does not en- 
tail these difficulties— although it is not entirely 
free from others. Here the change in magnetism 
of the crystal due to the introduction of hole elec- 
tron pairs injected by an external light source is 
measured directly. It is difficult to generate 
more than ~10*° carriers/cm*® by means of 
reasonable light intensities. Therefore, if one 
assumes the Landau-Peierls diamagnetism, one 
requires a detection sensitivity of at least 107" 
mass susceptibility units at room temperatures, 
which is ~100 times the sensitivity of the high- 
quality apparatus usually employed for meas- 
uring the y of solids. This sensitivity was in fact 
achieved in the present arrangement. The meth- 
od employed is basically the Faraday method, in 
which a small specimen is suspended in an in- 
homogeneous magnetic field. Our apparatus 
measures only the alternating component of the 
force produced when the sample is illuminated 
by light, modulated at the detection frequency of 
the apparatus (16 cps). The measurements are 

nevertheless “dc” since all pertinent processes 
are faster than the detection frequency. However, 


“de” drifts in the total x do not interfere with 
the observations. The detector is synchronous 
with the light modulation. 

The specimen used for all the measurements 
reported here was nearly intrinsic Ge (35 ohm- 
cm, n-type, at 300°K). The low-level minority . 
carrier mean life was 500 usec. The dimen- 
sions were 2X1x0.2 cm. For most of the meas- 
urements the specimen was etched so as to give 
low surface recombination velocity, s. A spe- 
cimen adjacent (in the parent crystal) to the one 
used for the x runs was used for observations of 
the photoconductivity. The photon flux and in- 
jected carrier density were determined by a 
combination of photoconductive decay lifetime 
measurements and steady-state photoconductivity. 

The dependence of the photoconductivity on 
light intensity was linear at room temperature, 
up to an absorbed photon flux of 3 x10'*/cm’ sec. 
At higher photon fluxes, a square root depend- 
ence set in, implying a (An)~! dependence of the 
mean life, 7. The dependence of the carrier 
density on the light level may easily be under- 
stood in terms of the Shockley- Read theory ap- 
plied to an intrinsic semiconductor.‘ 

Hole-electron pairs generated by light may 
give rise to two effects. The first is the diamag- 
netic susceptibility, x 7 of the added carriers; 
the other is a macroscopic loop current which 
produces a magnetic moment. Although loop 
currents could be eliminated through the use of 
pulverized specimens, none were employed be- 
cause of the difficulty in eliminating surface ef- 
fects and estimating the correct photocarrier 
density. 

To obtain the magnetic susceptibility, therefore, 
it is necessary to separate and remove that part 
of the force signal due to macroscopic moments 
from the total signal. The origin of the macro- 
scopic moments is the photoelectromagnetic 
(PEM) effect which may arise when there is a 
concentration gradient of carrier pairs in the 
presence of a magnetic field. The PEM short- 
circuit current, 7 s? is proportional to ¥4q*H, 
where Va is the diffusion current associated with 
the concentration gradient. When no electrical 
connections are made to the specimen, 7 ® flows 
within it. This circulating current gives rise to 
a magnetic moment, Eppyo ¥4*H, since the 
current does not cancel locally but only when 
summed over the profile of the entire specimen.°® 
The ipEM results in an observable force pro- 
portional to yu m’sradH. The effect has been 


predicted by Aigrain and Garreta® and implied by 
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Frenkel.’ ‘4 
When Va points into the sample, the force is 
paramagnetic in sign. When Va points from the 
interior outward, it is diamagnetic. The latter 
condition may be achieved either by treating the 
surfaces in such a way that the carrier recom- 
bination is higher there than in the sample in- 
terior, or by employing nonuniform illumination 
so that a component of Fe is perpendicular to H. 
The diamagnetic effect may occur when the pri- 
mary carrier diffusion is along H and thus does 
not produce a normal PEM effect. The PEM 
effects were exhaustively studied, both for their 
intrinsic value and so that a correct separation 
between them and the susceptibility could be 
performed. The short- circuit PEM current, i, 
which is proportional to CpEM was also meas- 
ured in the experimental geometry in order to 
permit comparison of its characteristics with 
the information from the force measurements. 
We found that if a specimen with etched, low 
surface-recombination-velocity surfaces is sus- 
pended in the inhomogeneous field and illumi- 
nated with very uniform light, two types of forces 
are seen. If \¥ 7H =vgH (90° orientation), a 
paramagnetic signal is observed, whereas when 
\¥7xHI =0 (0° orientation), the signal is diamag- 
netic. Both signals depend on the temperature, 
T, the paramagnetic signal rising and the dia- 
magnetic signal falling with increasing T. It was 
possible to resolve the 96° paramagnetic signal 
into two components having opposite signs. The 
paramagnetic component exhibited the character- 
istic PEM drop-off of quantum efficiency at long 
light wavelengths, as well as the same light in- 
tensity and (small) temperature dependence of 
the measured PEM ig As the incident light be- 
comes more penetrating at long wavelengths, the 
quantum efficiency of the normal PEM effect 
decreases sooner than that of the photoconductiv- 
ity due to the reduction of the carrier concentra- 
tion gradient. The second component of the force 
signal was diamagnetic and depended on T, light 
intensity, and wavelength in exactly the same 
way as the 0° effect. Thus the total general force 
signal was found to consist of one component due 
to the normal PEM effect, and a second one 
which was invariant to the relative orientation of 
incident light (or v4) with respect to H. Since a 
can be generated in the 0° orientation only 
whE™ nonuniform light and cannot remain invariant 
under a rotation of the sample through 90°, and 
since the observed T dependence of the diamag- 
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netic signal is much stronger than that of any 
PEM effect, we conclude that the diamagnetic 
signal originates in light-induced changes of the 
susceptibility, x,. In the 90° orientation the 
PEM and anmnine 

ly equal at 300°K and with a photon flux of 1.3 


x10'* photons/cm’? sec, whereas in the 0° orien- 
n 


tation only the susceptibility signal is observed, 
We have also shown, by means of a blackened 
specimen, that the diamagnetic signal is neither 
a manifestation of the Nernst effect, nor of syn- 
chronous temperature-induced fluctuations of 
the total x. 

The outstanding characteristic of y q 38 its 
magnitude; it is 400 times larger than the 
expected from summing the contributions of the 
photogenerated free carriers, assuming the 
Landau- Peierls theory and the cyclotron reso- 
nance effective masses averaged in the usual 
way.® The x, measured at 300°K for an incident 
photon flux J=1.3 x10'*/cm* sec, resulting in 
7x10"* pairs/cm’, is 1.02 x10~* mass suscepti- 
bility units. Of equally fundamental interest is 
the dependence of x, on J. Over the temperature 
range from 300°K to 420°K, x, is linear in/, 
whereas the light-generated pair concentration, 
4(an+ Ap), varies as J”? at 300°K and approaches 
linearity at the higher temperatures. In these 
experiments 2 x10'* <7 <10"* photons/cm’ sec. 
Furthermore, Xq drops quite rapidly as the tem- 
perature increases, although An + Ap does not. 
According to a preliminary empirical equation 
which fits all the data obtained so far, the drop- 
off in yj, for T in the vicinity of 400°K, is given 
by x, = const xIx (0.0025 ev/KT)(1 - 300/T)". 
This rather suggestive relation has not yet been 
tested at low T. The wavelength dependence of 
xX, is the same as that of the photoconductivity. 
It should be pointed out here that although the 
measured susceptibilities reported in this paper 
are anomalously large, they do not in fact di- 
rectly contradict any previous measurements. 
The conflict occurs only with the theoretical ex- 
trapolation of those measurements to the case of 
carrier pairs at our concentrations and tempera- 
tures. 

It is clear that xy, and An+ Ap are not linearly 
related. In the low-temperature range, where 
An+Ap>>mo+Po, the equilibrium carrier density, 
Xz» is apparently proportional to (An +p)’, 
whereas at the high T, where An+ Ap=no+Po: 
linear dependence is approached. It is therefore 
not unreasonable to suppose that the observed 
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anomalously large x , is due to some correlation 
or exciton effect, which would be expected to be 
proportional to a quadratic form in An, Ap, no, 
and Po: 

There exists the rather remote possibility 
that the light affects the large paramagnetic com- 
ponent of the lattice susceptibility directly. 
However, all processes which do not in some 
way relate to correlation or to another break- 
down in the usual theory may almost certainly 
be ruled out by the magnitude of the observed 
Le Further experiments, using extrinsic sam- 
ofes, together with measurements at lower tem- 
peratures, are expected to yield the answers to 
many of these questions. 

We wish to acknowledge the help of B. Tomp- 
kins in the construction of the apparatus, as 
well as several very helpful conversations with 
M. A. Lampert, J. Loferski, and A. Rose. 
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PIEZOELECTRIC PRODUCTION OF 
MICROWAVE PHONONS 


E. H. Jacobsen 
General Electric Research Laboratory 
Schenectady, New York 
(Received February 6, 1959) 


Recently Baranskii,' and B6mmel and Drans- 
fel have reported the production and propaga- 
tion of ultrasonic waves in crystalline quartz, 
excited by electromagnetic fields, in the fre- 
quency range of 10° to 2.5x10° cps. 

We have been able to extend the frequency, at 
which ultrasonic wave propagation is detectable, 
to the vicinity of 10000 Mc/sec. Specifically, 
we have successfully propagated longitudinal 





elastic waves at low temperatures in quartz bars, 
with optically flat ends and of various dimen- 
sions, placed between two re-entrant microwave 
cavities tuned to the transmitter frequency. 
Figure 1 shows an oscilloscope trace of the 
following experiment carried out in a liquid he- 
lium bath. A pulse of microwave energy at 9370 
Mc/sec is fed to one cavity, whereupon it is 
piezoelectrically converted to an elastic wave of 
the same frequency which travels along the 
quartz rod at the characteristic sonic velocity. 
Upon arrival of the wave at the receiver cavity, 
a part of its energy is reconverted to a radio- 
frequency signal and is detected and displayed 
on the oscilloscope, appearing as the first echo 
pulse. The remaining sound wave energy is re- 
flected at the receiving end of the rod, traveling 
back to the transmitting end whereupon it is 
again reflected to the receiver cavity to produce 
another echo pulse, and so on, until the ultra- 
sonic wave train dies out. It will be observed 
that, whereas the first echo pulse is detected at 
a time 7 after the direct transmitter pulse, the 
second and all subsequent pulses are separated 
by a transit time of 27 since these have made a 





FIG. 1. Oscilloscope trace of the receiver output 
registering pulses of microwave power (9370 Mc/sec) 
acoustically delayed by a quartz rod 3 cm in length and 
submerged in liquid helium at 4.2°K. The lower trace 
is identical with the one above except for an expanded 
time base. Time increases from left to right. Note 
that the first echo appears at a time 7 after the direct 
transmitter pulse (first pulse on photograph) whereas 
all subsequent echoes are spaced at intervals of 27 (see 
text). 
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complete round trip. This feature is seen more 
clearly in the lower trace of Fig. 1, which ex- 
hibits the pulses on an expanded time base. 

The dimensions of the quartz rod used in this 
experiment were 3 cm (length) by 0.3 cm (diam- 
eter), with the x axis along the length of the rod. 
The characteristic time (27) between pulses was 
about 11 microseconds, and the pulse duration, 
two microseconds. The peak rf power incident 
upon the transmitter cavity was 25 watts, de- 
rived from a 4J50 magnetron. 

In Fig. 1, the first dozen or so pulses appear 
to have nearly the same amplitude due to satura- 
tion of the receiver. After correction for this 
effect, the amplitude of each successive pulse 
is as shown in Fig. 2, where the third echo pulse 
is normalized to 100. The pulse decay envelope 
is seen to be more complicated in detail than 
the expected simple exponential drop-off. The 
details of the decay curve appear to be a function 
of the microwave electric field configuration at 
the point of coupling in the transmitting cavity 
and of the particular quartz rod employed. Thus 
the irregular decay may be caused by partial 
coupling to other elastic modes of the rod in 
such a way that energy is alternately exchanged 
between them. This detail is under further in- 
vestigation. 

When the temperature is varied from 1.8°K to 
20°K, it is found that the average attenuation 
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increases somewhat more rapidly than linearly 
with temperature. At 77°K no ultrasonic wave 
propagation was observable. It is expected that 
umklapp processes will become effective in 
attenuating the wave in the temperature range 
between 20°K and 77°K. 





'K. N. Baranskii, Doklady Akad. Nauk S.S.S.R. 114, 
517 (1957) [translation: Soviet Phys. Doklady 2, 237 
(1957)]. 7 

7H. E. Bommel and K. Dransfeld, Phys. Rev. Lett. 
1, 234 (1958). 





MAGNETIC OSCILLATIONS OF ULTRASONIC 
ATTENUATION IN A COPPER CRYSTAL 
AT LOW TEMPERATURES* 


R. W. Morse and J. D. Gavendat 
Department of Physics, 
Brown University, 
Providence, Rhode Island 
(Received February 24, 1959) 


In an earlier communication! we reported an 
apparently oscillatory dependence of ultrasonic 
attenuation in a polycrystalline sample of copper, 
this occurring at helium temperatures where the 
electron mean free path (1) was comparable to the 
ultrasonic wavelength (A). The results were ex- 
plained in terms of resonant conditions between 
the electron orbit diameter and the spatially per- 
iodic fields carried by the wave. It was shown 
that the average Fermi momentum calculated 
from such a picture was consistent with a value 
for copper of one free electron per atom. These 
measurements, however, were deficient in sev- 
eral obvious respects. Only one or two maxima 
and minima were discernible, and a polycrystal 
was used because at that time a single crystal of 
sufficiently long mean free path was not available. 
Furthermore, a subsequent theoretical analysis 
by Rodriguez” has suggested that such an oscil- 
latory effect would not be expected (at least from 
a classical Boltzmann equation analysis), and so 
it was proposed that the observed effect was per- 
haps a consequence of the polycrystalline nature 
of the sample. 

Here we report some more recent measure- 
ments made in a very pure copper single crystal 
which verify that a pronounced oscillatory effect 
is indeed found. Moreover, this effect is signifi- 
cantly anisotropic, suggesting that ultrasonic 
methods should be useful in obtaining information 





~ 











Vo 


abo 


lon; 


The 
but 
per: 
min 


puls 
of tl 
1.28 
with 
bett 


wav 
rect 
nor! 
rota 
fect 

plott 
1(a) 


+6- 


+2- 





(db/cm) 
4 ° 
i ae 


-44 


RELATIVE «<< 


-6-4 





“10 4+— 
is 


FIG, | 
tudinal | 
(a) is fo 
the [110 
to that f 
downwa: 








9 


a 


tal 
ct 
ifi- 


tion 








VoLuME 2, NUMBER 6 


PHYSICAL REVIEW 


LETTERS Marcu 15, 1959 








about details of the electronic properties of metals. 

Attenuation measurements were made of both 
longitudinal and shear waves at frequencies be- 
tween 15 and 75 Mc/sec and in magnetic fields of 
various orientations at strengths up to 6500 gauss. 
The temperature was varied between 1.0° and 4.2°K 
put the results were found to be essentially tem- 
perature -independent, indicating that / was deter - 
mined by impurity scattering. The attenuation 
was measured by observing a 1-ysec ultrasonic 
pulse reflected between two highly parallel faces 
of the sample, these being separated by about 
1.25cm. The faces of the samples were oriented 
with respect to the crystalline axes to an accuracy 
better than 3°. 

Figure 1 shows the results for a longitudinal 
wave of 75 Mc/sec propagated along the [001] di- 
rection. The magnetic field was applied in a plane 
normal to the propagation direction and could be 
rotated about this axis. Since the oscillatory ef- 
fect should be periodic in H™, the results are 
plotted in Fig. 1 on a reciprocal scale. Figure 
1(a) shows data for H along either a [100] or a 
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FIG, 1. Attenuation along [001] of 75-Mc/sec longi- 
tudinal wave vs AH plotted on reciprocal scale. Curve 
(a) is for H in the [100] direction, and (b) is for H in 
the [110] direction. Attenuation is measured relative 
'o that for H=0, with the (b) curve arbitrarily displaced 
downward by 6 db/em. 


[010] direction. About seven maxima and minima 
are found which are quite accurately periodic in 
H™, Only slightly different results are found for 
H within about + 15° of these directions. However, 
with H in the vicinity of the 45° directions, such 
as [110], the oscillations are changed consider - 
ably, a typical result being shown in Fig. 1(b). 

Although no satisfactory theory for such oscil- 
latory behavior now exists, it is of interest to ob- 
serve that the results are consistent with a model 
of the Fermi surface proposed by Pippard based 
on his observations of the anomalous skin effect. 
In the ultrasonic effect one would expect the great- 
est contribution to come from those electrons in 
an extremal cross section normal to the direction 
of H. Pippard suggests that the Fermi surface 
in copper touches the zone boundary in the [111] 
directions, but is nearly spherical elsewhere. 
Thus for H along the principal axes the relevant 
electrons would lie on a nearly circular cross sec- 
tion of the Fermi surface. This is the orientation 
where the well-defined periodic effect is found. 
On the other hand, when H is along a 45° direc- 
tion the relevant electrons are on a cross section 
that is broken up into four segments by the zone 
boundaries. Experimentally, it is in these direc- 
tions that one sees a more complicated effect. In 
particular, another oscillation of a different per- 
iod seems to enter. A maximum of this other 
period is designated by x in Fig. 1(b). 

One can make an estimate of the Fermi momen- 
tum (p,-) from the measurements shown in Fig. 1. 
If it is assumed that between adjacent maxima the 
average orbit radius changes by A/2, and that 
Y=p,c/eH, where p, is the momentum component 
perpendicular to H, then one finds p, =1.16x107® 
g cm/sec for H along a principal axis. It seems 
likely that it is some average of p , that matters 
and not the extremal value p,. If the assumption 
of a spherical shell is used, then Pp, =(1/4)p¢, or 
by =1.47x 10°*° g cm/sec. According to Pippard’s 
estimate, p, in the extremal cross section per- 
pendicular fo [100] varies between 1.39 x 107% 
and 1.44107 g cm/sec. 

Contrary to the results in the polycrystal, the 
attenuation does not seem to approach zero as 
H becomes large (at least in all orientations). A 
similar effect was noted earlier in a tin single 
crystal.* The attenuation at high values of the field 
(i.e., for 7 <A/2) also is found to be highly aniso- 
tropic in copper, this dependence being shown in 
Fig. 2. The angular variation shows a fourfold 
symmetry which is distorted somewhat, probably 
because of unavoidable misalignment with respect 
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FIG, 2. 
tudinal wave vs the direction of H about the [001] axis, 
where H=6500 gauss. The direction (a) corresponds 
to [100] and (b) to [110]; these are the directions in 
which the curves shown in Fig. 1 were taken. The 
attenuation scale is relative to that at H=0. 


to the magnet. The interpretation one should place 


on this polar pattern is not clear. Rodriguez’s 


theory predicts in this situation that the high-field 


attenuation should become proportional to w?7 


(where w is the angular frequency of the wave and 


7 is the electron relaxation time). Thus the pat- 
tern shown in Fig. 2 may reflect anisotropy of 7 
as well as of the Fermi surface. 


Preliminary measurements also have been made 
of the attenuation of a plane polarized shear wave 


with H along the direction of propagation. As 
pointed out by Kjeldaas, this configuration is in- 
teresting because certain groups of electrons, 
selected by the value of the magnetic field, re- 
main in resonance with the periodic fields.*> The 


results at the lower frequencies agree qualitative- 


ly with Kjeldaas’ predictions except that the ex- 


perimental attenuation initially increases with H. 

The attenuation at high frequencies does not dimin- 
ish with increasing H as rapidly as the theory pre- 
dicts. However, this may be due to a magnetically 


induced rotation of the plane of polarization, an 
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effect predicted by Kjeldaas. These results, as 
well as measurements in other crystallographic 
directions, will be discussed at more length in 
a subsequent paper. 





“Partially supported by a contract with the U. S. Air 
Force through the Office of Scientific Research of the 
Air Research and Development Command. 

TEdgar Lewis Marston Fellow. 

‘Morse, Bohm, and Gavenda, Phys. Rev. 109, 1394 
(1958). 

*§. Rodriguez, Phys. Rev. 112, 80 (1958). 

34, B. Pippard, Phil. Trans. Roy. Soc. (London) 
A250, 325 (1957). 

‘Morse, Bohm, and Gavenda, Bull. Am. Phys. Soc, 
Ser. Il, 3, 44 (1958). 

'T, Kjeldaas, Jr., Bull. Am. Phys. Soc. Ser. I, 
3, 180 (1958); also, Westinghouse Scientific Paper 
6-94439-1-P8, 1958 (unpublished). 





SOME EFFECTS OF STRAIN IN THIN 
SPECIMENS ON ABSORPTION 
MEASUREMENTS AT LOW TEMPERATURE 


G. G. Macfarlane, T. P. McLean, 
J. E. Quarrington, and V. Roberts 
Royal Radar Establishment, 
Great Malvern, England 
(Received February 16, 1959) 


Direct transition excitons’’* and the magneto- 
optical effect?»* have recently been investigated 
in Ge by making transmission measurements on 
polished plates 10 microns and less in thickness, 
over the range from room temperature to 1.5°K 
under high spectral resolving power. We have 
already drawn attention’ to discrepancies be- 
tween our measurements and those of the Lin- 
coln group,* which are much greater than the 
quoted limits of error. What we believe to be 
the explanation of these differences is now pre- 
sented. It is emphasized that in order to get re- 
sults which are accurately characteristic of the 
specimen, it must be mounted in a strain-free 
manner. 

The discrepancies are twofold: 

(a) At low temperatures the absorption line 
due to excitation of electrons into the ground 
state of the direct transition exciton is consist- 
ently observed by the Lincoln group at a higher 
energy than by us. 

(b) The direct transition energy gaps E, de- 
duced from the observations, which are in ex- 
cellent agreement at room temperature, ‘ diverg¢ 
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progressively on cooling, the Lincoln values 
being the higher. 

Since (b) involves differences of interpretation, 
we focus attention on (a) which is entirely a 
matter of observation. 

With the willing cooperation of the Lincoln 
group we have been able to make measurements 
on one of their specimens. For convenience this 
was done at 77°K, and confirmed the position of 
the exciton line in the Lincoln group’s specimen 
to be in good agreement with their own measure- 
ments? (within ~0.001 ev). We therefore discount 
errors on either part in spectrometer calibra- 
tion and difference in specimen temperature, and 
turn attention to the specimens themselves. The 
difference between them may be intrinsic or due 
to their experimental environment, but as we 
have noticed no significant variation between 
specimens taken from different crystals, the 
latter is the more likely. 

A probable explanation was in terms of strains 
inthe specimens, which would distort the band 
structure and show up as a shift or broadening 
of the absorption edge. Since the Lincoln speci- 
mens were glued to a glass backing they would 
certainly be under considerable strain on cool- 
ing, owing to differential contraction, although 
probably comparatively free from strain at 
room temperature, whereas ours being unbacked 
and freely suspended are strain-free at all tem- 
peratures. This would account for the progres- 
sive disparity between the observations on cool- 
ing. Furthermore our measurements and the 
Lincoln group’s are in fairly good agreement in 
the indirect transition region.”»> Here the ab- 
sorption coefficient is about 100 cm™ or less 
and thick specimens are used, which would be 
much more difficult to strain, even if backed. 

These ideas were tested as follows: 

(i) The Lincoln technique was duplicated as 
far as possible, i.e., a specimen (RRE1) glued 
to glass with Araldite® was prepared and meas- 
ured at 77°K. In this, the exciton line appeared 
ata higher energy in rough agreement with the 
Lincoln value. Perfect agreement would not be 
expected as the strain taken up by the specimen 
depends on a number of details such as the ex- 
pansion coefficient of the glass, the adhesive 
and its thickness, etc. It may be noted that 
wing crystal quartz instead of glass produced a 
considerably larger shift to higher energy. 

(ii) The position of the exciton line was ob- 
served in a freely mounted specimen (RRE2) 
vhich was then glued to a glass substrate and 


remeasured. The line was shifted to higher en- 
ergy as expected. These results are plotted in 
Fig. 1. In addition to the shift it is seen that 
backing the specimen also distorts the form of 
the transmission curve. In particular an extra 
band (marked a in Fig. 1) appears. This band 
was noted by the Lincoln group® but its presence 
was unexplained. Our measurements indicate 
that it is produced by strain in the specimens 
and shifts along the energy scale with tempera- 
ture in the same way as the absorption edge. It 
may well be linked with distortion of the band 
structure. 

We have already shown’ that the absorption 
coefficient of a freely mounted specimen follows 
Elliott’s’ theoretical predictions at energies 
above E,. This is not so for the backed speci- 
mens. 

(iii) Since the direction of shift was consistent 
with compressional strain, * a system was 
sought which would impose a tension on the spe- 
cimen when cooled, and should therefore produce 
a shift in the opposite direction. Assuming the 
backing to be the important factor, then this 
should be achieved by mounting a sample on 
fused silica.® This was tried and a shift of about 
0.005 ev occurred on cooling, in the direction of 
lower energy as expected. It is interesting that 
Rose-Innes® used a very similar system deliber- 
ately to strain samples of Ge and Si. He esti- 
mated that the stresses involved correspond to 
about 10° dyne cm™?. Taking Paul and Warsch- 
auer’s estimate® of the pressure coefficient of 
E, for Ge to be 1.2x10™“' ev/dyne cm™ we might 
expect a shift of around 1.2 x10™ ev for this 
specimen configuration if all the strain is taken 
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FIG. 1. Transmission ratio of various Ge specimen 
configurations at 77°K. 
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up by the specimen. This is about twice that 
actually observed. 

The results clearly show that thin specimens 
must be mounted in a strain-free fashion for 
this type of measurement if data characteristic 
of the unstrained material are desired. It fol- 
lows that cooling the specimen by conduction 
through a solid backing is unsatisfactory since 
it is practically impossible to arrange for the 
relative expansion to be zero. Recourse must 
be had to direct immersion in the cooling liquid, 
which brings its own difficulties, '® or immersion 
in thermal exchange gas.*'»** The latter method 
is used in this laboratory. These considerations 
will apply, of course, to other materials than 
Ge. 

Strain in the specimens will also affect mag- 
neto-absorption measurements. Clearly, a uni- 
form strain will shift the band edges and the 
adjacent Landau levels to roughly the same ex- 
tent. The magneto-absorption bands will thus 
experience an appreciable shift (of the same or- 
der as that of the exciton absorption), but it is 
unlikely that changes in their relative positions 
will be detectable by the spectroscopic means 
currently in use. The practical situation with a 
backed specimen is almost certainly one of non- 
uniform strain and will be more complicated. 

Finally it is emphasized that the condition of 
the specimens should be kept very much in mind 
when correlating data obtained from different 
experiments. 

We are indebted to the Lincoln group for lend- 
ing us two of their specimens and for communi- 
cating their results to us before publication. 
Our thin specimens were skillfully prepared by 
S. C. J. Brooke. This Letter is published by 
permission of H. B. M. Stationery Office. 





‘Macfarlane, McLean, Quarrington, and Roberts, 
Proc. Phys. Soc. (London) 71, 863 (1958). 

2Zwerdling, Lax, Roth, and Button, Phys. Rev. (to 
be published). In this paper the transmission ratio of 
a sample at 77°K is given and can be directly compared 
with out results. 

’Zwerdling, Lax, and Roth, Phys. Rev. 108, 1402 
(1957). 

‘In reference 1 we state that our value for Ey at 291°K 
should be regarded as less reliable than those at the 
lower temperatures. This is so, but should not be 
taken to imply that it is unreliable in the present con- 
text. The discrepancies which we are now discussing 
(~0. 006 ev) are more than ten times as great as the 
error in Ey. 

‘Macfarlane, McLean, Quarrington, and Roberts, 
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Phys. Rev. 108, 1377 (1957). 

A thermosetting adhesive similar to Allymer cement, 
We understand Allymer cement is used by the Lincoln 
group. 

"R. J. Elliott, Phys. Rev. 108, 1384 (1957). 

*w. Paul and D. M. Warschauer, J. Phys. Chem. 
Solids 5, 89 (1958). 

*A. C. Rose-Innes, Proc. Phys. Soc. (London) 72, 
514 (1958). ve 

1%V. Roberts, Colloquium Spectroscopicum Internatio- 
nale VI (Pergamon Press Ltd., London, 1956), p. 508, 
‘' McMahon, Hainer, London, and King, J. Opt. Soc. 

Am. 39, 786 (1949). 
vy. Roberts, J. Sci. Instr. 32, 294 (1955). 








MAGNETIC FORM FACTOR OF COBALT* 


R. Nathans 
The Pennsylvania State University, 
University Park, Pennsylvania, and 
Brookhaven National Laboratory, 


Upton, New York 
and 


A. Paoletti 
Comitato Nazionale Ricerche Nucleari, 
Rome, Italy, and 
Brookhaven National Laboratory, 
Upton, New York 
(Received February 26, 1959) 


We wish to report preliminary results on the 
magnetic form factor of the face-centered cubic 
phase of metallic cobalt, which indicate a devia- 
tion from spherical symmetry for its 3d elec- 
trons. These measurements represent a con- 
tinuation of studies on the angular dependence of 
the magnetic scattering in the transition elements 
using a polarized neutron beam spectrometer’; 
the results on iron and nickel have been already 
reported in reference 1. These studies were 
undertaken because of the current interest in the 
3d electron configurations in the transition ele- 
ments and the pertinence of magnetic form fac- 
tor data to such discussions. 

The sensitivity of the polarized beam technique 
in establishing the absolute magnitude of the 
magnetic scattering depends on the ratio of the 
nuclear scattering amplitude to the magnetic 
scattering amplitude at the different scattering 
angles. For cobalt, this ratio is very favorable 
and it is possible to obtain sizable differences in 
the reflectivities for beams of different polariza- 
tion even for large scattering angles, where the 
magnetic scattering is greatly reduced. For this 
reason every effort was made to reduce the in- 
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trinsic experimental uncertainties present in this 
type of measurement to a minimum so as to as- 
certain whether in cobalt the angular dependence 
of the magnetic scattering gives any definite in- 
dication of structure in the distribution of un- 
paired electrons. 

The method used to obtain the absolute magnetic 
scattering amplitude has been described in detail 
inreference 1. Briefly, the ratio of reflected 
intensities from a magnetized single crystal is 
measured for incoming neutrons polarized first 
parallel to the field on the crystal, and then anti- 
parallel. This polarization ratio for a particular 
Bragg reflection depends on the ratio p, pi’ ppp 
where ppp; and b,,, are the magnetic and nuclear 
structure factors for the (hk/) reflection. Since 
bye] is known, the measurement yields Pret: It 
is necessary to correct the observed values of 
the polarization ratio for the presence of extinc- 
tion in the test crystal and for the depolarization 
of the neutron beam as it traverses the crystal. 
The latter was measured by a method given in 
reference 1, and was shown to be less than 1%. 
Extinction, although small, was estimated by 
performing the measurements for crystals of 
different thicknesses and at two different neutron 
wavelengths (0.89A and 1.04A). 

The test crystal for the fec measurements con- 
sisted of a pillar crystal of the alloy Co, ..F ep. og, 
the presence of the iron serving to stabilize the 
cobalt in the fcc phase. All measurements were 
made at room temperature. The results of the 
measurements taken to date are shown in Table 
land plotted in Fig. 1. In calculating the mag- 
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FIG. 1. Magnetic form factor of face-centered cubic 
cobalt, 


Table I. Magnetic form factor of face-centered cubic 
cobalt. 








Reflection sin 6/A Pppi/ nei f 
(111) 0.244 1.21 +0.06 0. 781+0.040 
(200) 0.283 1.01 +0.06 0.652+0. 040 
(220) 0.397 0.57520. 005 0.37120. 003 
(311) 0. 467 0.41840. 003 0.27020. 002 
(222) 0.487 0.350+0. 002 0. 26620. 002 
(400) 0.563 0.286+0. 006 0.18420. 004 
(331) 0.613 0. 16620. 002 0.107+0.002 
(422) 0.689 0.09120. 006 0.059+0. 004 
(333) 0.731 0.05420. 006 0.03520. 004 
(440) 0.796 0.03440. 003 0. 022+0. 002 





netic scattering amplitudes we have used the 
values b= 0.25 x10" cm? and bye = 0-96 x 10°” 
cm. 

Of particular interest in the data is the failure 
of the form factor values for some of the higher 
(hkl) reflections to fall on a smooth curve, in 
particular the (311) and (400) reflections. In 
order to emphasize that these “bumps” are real, 
we have plotted in Fig. 2 the outer portion of the 
form factor curve on an expanded scale along 
with the estimated uncertainties in the measure- 
ments. These “bumps” have led us to suggest a 
deviation from spherical symmetry for the 3d 
electrons. As far as we know, these measure- 
ments give the first experimental indication of 
such a departure. A somewhat similar behavior 
to the above had also been noticed in the fcc 
nickel form factor. However, because of the low 
magnetic scattering amplitude of nickel, espe- 
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FIG. 2. Outer portion of the magnetic form factor 
curve of fec Co. 
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cially at the higher scattering angles, it was 
difficult to place these observed “bumps” on a 
more definite basis. 

Because of the change in the symmetry, it 
would also be of interest to compare the magnetic 
form factor of hexagonal close-packed Co with 
that of face-centered cubic Co. Hence we are 
continuing these measurements in order to refine 
the fcc data further and to determine the form 
factor of the hexagonal phase. A detailed report 
will be submitted upon the completion of this 
work. 

We are indebted to Dr. T. A. Kaplan and Dr. 

S. J. Pickart for profitable discussions of these 
results. We also wish to thank Dr. H. J. Williams 
and Dr. R. M. Bozorth of Bell Telephone Labora- 
tories for supplying us with the crystals. 





‘This work was performed under the auspices of the 
U. S. Atomic Energy Commission and the National 
Security Agency. 

‘Nathans, Shull, Shirane, and Andresen, J. Phys. 
Chem. Solids (to be published). 

2w. L. Roth, Phys. Rev. 110, 1333 (1958). 





LATTICE VIBRATIONS IN SILICON 
AND GERMANIUM 


B. N. Brockhouse 
Physics Division, 

Atomic Energy of Canada Limited; 
Chalk River, Ontario, Canada 
(Received February 13, 1959; revised 
manuscript received March 2, 1959) 


The dispersion curves of lattice waves travel- 
ling in the [001] directions of a silicon single 
crystal have been measured using neutron scat- 
tering. The frequency (v) and wave vector @) of 
each phonon, and the character of each branch, 
were determined as previously discussed.’ The 
complete results for the longitudinal acoustical 
(LA), longitudinal optical (LO), and the degener- 
ate transverse acoustical (TA) and transverse 
optical (TO) branches are shown in Fig. 1. 

Two series of experiments were performed 
using a new multiple axis crystal spectrometer 
at the NRU reactor. In the first, neutrons of 
1.900A wavelength were scattered from a single 
crystal of silicon at angles in the vicinity of 90° 
with the crystal oriented so that the wave vector 
of the outgoing neutrons (K’) lay very nearly 
along the line (2, 2, 3) to (2, 2, 6) of the silicon 
reciprocal lattice. The wave vectors of all pho- 
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tice vibrations propagating in the [001] direction of the 
reduced zone in silicon at room temperature. The 
heavy lines have the slope of the appropriate velocity 
of sound as found by ultrasonic measurements [H. J. 
McSkimmin, J. Appl. Phys. 24, 988 (1953)]. The 
vertical arrow indicates the presumed position (q) of 
the minimum of the conduction band (reference 13). 
The horizontal arrows indicate the frequencies of 
phonons at q, found from the indirect transition (ref- 
erences 10-12). 


nons observed in these experiments thus lay in 
the [001] direction of the reduced zone. In the 
second series, neutrons of 1.900, 1.350, and 
1.693 angstroms were used in the method of 
successive approximations,‘ in which phonons 
were measured in regions of reciprocal space in 
which their polarization, and optical or acousti- 
cal character, could be unambiguously deter- 
mined. All measurements were made at room 
temperature. The points are believed to define 
the level of all branches to about 3%. Values of 
important frequencies (units 10’ cps) are: Ra- 
man (q=0) 15.3+0.3, TO[001] 14.2+0.3, L[001] 
11.94+0.5, TA[001] 4.3540.15. 
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The measurements are in disagreement with 
x-ray measurements of Learn.* The Raman fre- 
quency is in agreement with a prediction* 15.45 
x10 cps based on the far infrared spectrum, 
and with an independent neutron measurement of 
15.210" cps by Palevsky, Hughes, Kley and 
Tunkelo.® The dispersion curves are in agree- 
ment with predictions*»*® from specific heat, 
elastic constant, and far infrared measurements, 
that the frequencies for silicon are similar to 
those for germanium at similar positions in the 
zone but higher by a factor of about 1.75. Six 
figures of merit obtained from these data gave* 
an average ratio 1.745 with a standard deviation 
of 0.04. The average ratio of the silicon fre- 
quencies reported herein to the frequencies of 
germanium measured by Brockhouse and Iyengar’ 
at the same positions in the reduced zone is 1.73 
+0.04. If exact homology held between silicon 
and germanium, the harmonic force constants 
would be proportional to the inverse square of 
the lattice constant (a), and the frequencies 
therefore proportional to (Ma*)~””, where M is 
the mass of the atom. Thus the frequency ratio 
would be 1.675, slightly but significantly lower 
than the measured ratios. The expected fre- 
quency ratio can also be computed from the 
melting temperatures using the Lindemann law,’ 
leading to a value of 1.96. The lack of complete 
homology thus seems to be connected with an 
abnormally low melting point in germanium 
(compared with silicon). 

It follows from the approximate homology that 
in silicon, as in germanium, either long-range 
forces between atoms are important or the 
Born-Von Kdrman theory is inadequate.'» *: °. 

The [001] direction was selected for initial 
study because the minimum of the conduction 
band lies in this direction, and it was hoped that 
the wave vector (q,) of the minimum could be 
located by comparison of these results with other 
types of measurements which are dependent on 
the value of G5. Measurements of the infrared 
absorption spectrum by Macfarlane, McLean, 
Quarrington, and Roberts’° have been interpreted 
to give values for the frequencies at q, of 4.4, 
14.0, 21.9, and 29.5 x10" cps. Measurements of 
the spectrum of recombination radiation by 
Haynes, Lax, and Flood" originally gave fre- 
quencies a little lower than those obtained by 
Macfarlane et al., but more recent measure- 
ments are in agreement.” The two higher fre- 
quencies of each set are clearly inconsistent 
with the results of Fig. 1 but could arise from 


double phonon processes as suggested‘ earlier. 
From measurements of electron nuclear double 
resonance, Feher’* has found that q)/q,,45 = 9-85 
+0.03. If this value is accepted, then the two 
lower frequencies of each set must be TA and 
probably TO. The LO and LA transitions are 
presumably weak and overshadowed by the very 
strong TO transition. 

The measurements of lattice frequencies of 
germanium of Brockhouse and Iyengar,’ dis- 
cussed above, and those of Ghose, Palevsky, 
Hughes, Pelah, and Eisenhauer“ are in disagree- 
ment for the TA and especially for the longitudi- 
nal branches near the zone boundary in the [001] 
direction. For L[001] Brockhouse and Iyengar 
quoted (6.9+0.4) x10"? cps, while the correspond- 
ing frequency derived from the measurements of 
Ghose et al. is (5.3+0.4) x10" cps. In an effort 
to resolve this question we have remeasured’® 
the [001] longitudinal branch in germanium using 
neutrons of longer wavelength (1.917A) than those 
used by Brockhouse and Iyengar (1.52A), and 
somewhat higher resolution. The two types of 
experiments discussed above for silicon were 
performed, and gave results in agreement with 
one another. The dispersion curve (nine points) 
obtained agreed with that of Brockhouse and 
Iyengar, the frequency of the L[001] phonon being 
(7.0+0.3) x10" cps. 

The author wishes to thank Dr. J. R. Haynes, 
Dr. M. Lax, Dr. J. Burton, and Dr. M. Tanen- 
baum of Bell Telephone Laboratories for the 
loan of the crystals. 





'B. N. Brockhouse and P. K. Iyengar, Phys. Rev. 
111, 747 (1958). 

*See Fig. 4 of reference 1. 

3A, Learn, thesis, Massachusetts Institute of Tech- 
nology, 1958 (unpublished). See also H. Cole and E. 
Kineke, Phys. Rev. Lett. 1, 360 (1958). 

‘B. N. Brockhouse, paper given at the International 
Conference on Semiconductors, Rochester, New York, 
August 18-22, 1958 [J. Phys. Chem. Solids (to be pub- 
lished)]. 

*Palevsky, Hughes, Kley, and Tunkelo, this issue 
(Phys. Rev. Lett. 2, 258 (1959)]. 

‘Similar homology arguments have been presented 
in the past; see references quoted in reference 4 and 
F. Herman (private communication). 

"See N. F. Mott and H. Jones, Properties of Metals 
and Alloys (Oxford University Press, London, 1936), 
p. 13. 

8¥, Herman, paper given at Conference on Semi- 
conductors (see reference 4). 

°M. Lax, Phys. Rev. Lett. 1, 133 (1958). The argu- 
ments of Lax on silicon are still applicable even though 
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the actual phonon frequencies are greatly different 
from those he used, which were those given by the in- 
frared workers. 

10Macfarlane, McLean, Quarrington, and Roberts, 
Phys. Rev. 111, 1245 (1958); paper given at Conference 
on Semiconductors (see reference 4). 

‘tHaynes, Lax and Flood, Bull. Amer. Phys. Soc. 
Ser. I, 3, 30 (1958); paper given at Conference on 
Semiconductors (see reference 4). 

123. R. Haynes (private communication). 

8G. Feher, paper given at Conference on Semicon- 
ductors (see reference 4). 

“Ghose, Palevsky, Hughes, Pelah, and Eisenhauer, 
Phys. Rev. 113, 49 (1959). 

These experiments were performed after the initial 
submittal of this Letter. 





LATTICE VIBRATIONS IN SILICON BY 
SCATTERING OF COLD NEUTRONS* 


H. Palevsky, D. J. Hughes, W. Kley,t 
and E. Tunkelo 
Brookhaven National Laboratory, 
Upton, New York 
(Received February 18, 1959) 


The use of subthermal or “cold” neutrons for 
investigating the lattice vibrations in solids has 
proved to be of increasing value in the past 
several years. Upon absorption of a phonon (the 
quantum of lattice vibration energy) while scat- 
tering, the neutron energy changes by such a 
large amount that the phonon energy can usually 
be accurately fixed. The method is particularly 
useful for the optical vibrations, as shown by 
the first measurements’ of high-frequency modes 
in zirconium hydride and germanium, cases for 
which infrared work had given inconclusive re- 
sults. In the light of similar uncertainty con- 
cerning the optical (Raman) frequency in silicon, 
we have determined it, as well as one acoustical 
branch, by use of cold-neutron scattering. 

Because of the great interest in the electronic 
properties of silicon, a number of infrared meas- 
urements have been performed, which have been 
used to infer the Raman frequency, the optical 
vibration of zero phonon momentum. For ex- 
ample, Collins and Fan,” on the basis of infra- 
red absorption measurements, conclude that hv 
for the Raman line is 0.076 ev. Haynes, Lax, 
and Flood,* from a measurement of recombina- 
tion lines from infrared excitation, obtain optical 
lines between 0.08 and 0.12 ev. Macfarlane, 
McLean, Quarrington, and Roberts,‘ on the basis 
of infrared absorption edge measurements, re- 
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port similar values, 0.09 and 0.12 ev. Recently, 
Cole and Kineke® have reported that the unpub- 
lished x-ray results of A. Learn give a value of 
0.071 ev. On the theoretical side, Hsieh, ® using 
the theory of Smith,’ has calculated the optical 
vibration spectrum from the measured force 
constants and obtained a Raman value of 0.077 
ev, using only first neighbor interactions, anda 


value of 0.10 ev when the effects of second neigh- 


bors are included. These discordant results in- 
dicate the need for an unambiguous value of the 
optical line, such as that obtainable from cold- 
neutron scattering. 

In the cold-neutron technique as developed at 
Brookhaven, ® the neutron beam incident on the 
scatterer is produced by filtration through 8 
inches of beryllium. No Bragg (elastic) scatter- 
ing can occur for the low-energy (~ 0.005 ev) in- 
cident neutrons, and those scattered at 90° (the 
position of the detector) are of necessity those 
that have gained energy by interacting with the 
lattice vibrations. Since the energy of the optical 
vibrations of silicon is in the neighborhood of 
0.05 ev, the neutrons experience a tenfold in- 
crease in energy when scattered (Fig. 1). Itis 
this large change in energy, giving directly the 
phonon energy, that is a distinct advantage of 
the cold-neutron method. 

In the present measurements, the incident neu- 
tron beam was perpendicular to the (100) plane 
of the silicon single crystal. The energy of neu- 
trons scattered at 90° was determined by meas- 
uring their time of flight over a 5-meter flight 
path, utilizing a slow chopper. Figure 1 shows 
a typical set of experimental data, plotted in 
terms of time of flight. The sharp peak at high 
energy represents the absorption of one quantum 
of vibrational energy from an optical mode, and 
the one at lower energy corresponds to an acous- 
tical mode. Data similar to Fig. 1 were obtained 
for a number of orientations of the crystal as it 
was rotated about the [100] axis. The results 
give energies for phonons of various momenta q, 
with q limited to the (100) plane of reciprocal 
lattice space. Points for q in the [110] direction, 
extrapolated to q=0, give the energy of the Ra- 
man line as 0.063 + 0.003 ev (phonon frequency, 
v, of 15.2x10" cps); and at the edge of the Brill- 
ouin zone (maximum 4Q) in the [110] direction the 
transverse acoustic branch has an energy of 
0.020+ 0.002 ev (v=4.84x10" cps). Thus our 
value for the Raman line is not in agreement with 
the infrared measurements nor with the Hsieh 
calculations. It is possible, however, to reinter- 
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FIG. 1. 


Counting rate as a function of time of flight for neutrons scattered at 


90° by Si. The arrow indicates the position of the incident neutrons on the same 
scale. The peaks marked Ey and Ey, are the energies of scattered neutrons 
which have gained one quantum of vibrational energy from the optical and acous- 


tical modes, respectively. 


pret the high-energy infrared results in terms of 
combination lines so that they would be consistent 
with the neutron measurements. 

Our purpose in measuring one acoustical branch, 
in addition to the Raman line, was to check the 
validity of a simple model suggested about two 
years ago by Herman, ° who argued that in the 
first approximation the vibration frequencies of 
Si and Ge should be inversely proportional to the 
square root of the masses of the vibrating atoms, 
for they both have the same valence electron 
structure. This ratio is 1.61, and a comparison 
of the present results with our Ge measurements’? 
gives a ratio of 1.69+0.10 for the Raman frequen- 
cies @=0) and a ratio of 1.61+0.15 for the trans- 
verse acoustic phonons at the edge of the Brill- 
ouin zone in the [110] direction. Thus an ex- 
tremely simple physical model serves to corre- 
late the frequencies of Si and Ge within the pre- 
sent experimental error, for both optical and 
acoustical modes. The failure of Hsieh’s calcu- 
lation of the vibration frequencies, and the pro- 
portionality of Si to Ge frequencies, indicate 
that for Si, as well as Ge, there are long-range 
forces present. 


Note added in proof.—The results for the Raman 





frequency in this Letter and those of Brockhouse 
in the accompanying Letter (using neutrons of 
higher energy) are in excellent agreement, for 
both silicon and germanium. There are unex- 
plained differences for the acoustic branches, 


for which some of our frequencies are higher, 
and some lower than his. The differences, how- 
ever, do not affect the conclusion that the ordi- 
nary Born—Von Karman theory using short- 
range forces only is inadequate, which follows 
from either set of results. 





*Work performed under the auspices of U. S. Atomic 
Energy Commission. 

TGuest scientist from University of Freiburg, Frei- 
burg, Germany. 

tGuest scientist from Finland Institute of Technology, 
Helsinki, Finland. 
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CONTRIBUTION OF THE NUCLEUS TO THE 
SPECIFIC HEAT OF RHENIUM* 


P. H. Keesom and C. A. Bryant 
Purdue University, Lafayette, Indiana 
(Received February 25, 1959) 


Measurements of the heat capacity of a Re in- 
got? between 0.37 and 4.2°K show an anomalous 
contribution to the specific heat which is pro- 
portional to T~?. This temperature dependence 
suggests that the contribution is the high-tem- 
perature part of a Schottky-type anomaly, arising 
from excitations among a set of discrete energy 
levels. Figure 1 shows the ratio of the total 
specific heat to the absolute temperature, C/T, 
plotted as a function of T. The specific heat of 
the Re measured in zero magnetic field can be 
expressed as 


C,=yT+aT*+AoT”’, T>1.85°K (1) 
C,=C,,+aT%+AgT?, T<1°K (2) 


where yT is the electronic specific heat in the 
normal state, aT® is the specific heat of the 
lattice, and A,T~* is the anomalous contribution. 
Cos is the electronic specific heat which is ex- 
pected, according to the Bardeen, Cooper, and 
Schrieffer theory of superconductivity,? to be 


given by 


Cog =7T,aexp(-bT,/T), T/T >2. (3) 


Because aT® is small throughout and A,7~ is 
small at the higher temperatures, the constants 
can readily be adjusted by successive approxi- 
mations to fit the data. Accepting the validity of 
Eq. (3), we find A,=0.052 40.002 millijoule 
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FIG. 1. Ratio of the specific heat of rhenium to ab- 
solute temperature as a function of temperature in 
zero magnetic field (e) and in 1000 oe (4). The solid 
line is C,./T vs T. 
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deg/mole. The other constants are y=2.31+0.02 
millijoule/mole deg” and a =0.026+0.002 milli- 
joule/mole deg*, from which it follows that the 
Debye temperature 6,=417+10°K. With this 
choice of constants, the data appear to be con- 
sistent with the second law of thermodynamics, 
i.e., with 

z. T. 


c c 
a‘ o2 e 
(/ Cog? ar | vat) 
) fr) 


within the limits of experimental accuracy. 
Setting T,,=1.70°K in Eq. (3), we obtain a =8.5 
and 5 =1.46. 

To quench the superconductivity and observe 
the magnetic field dependence of the anomaly, we 
remeasured in a field of about 1000 oersted. 
The resulting specific heat is also plotted in Fig. 
1, and can be expressed as yT+a@7T°+A,T” 
with A,,=0.06+0.01 millijoule deg/mole (because 
of fluctuations in the applied field, this measure- 
ment was less accurate). Ay does not appear 
significantly different from A,, and we conclude 
that the anomaly is influenced by the applied 
field less than should be expected if it were 
caused by repopulation of electronic energy 
levels. For, since no deviation from the T” de- 
pendence of the anomaly is observed down to 
0.37°K, the characteristic temperature of the 
zero-field splitting is certainly less than 0.37 
deg. Calculations for a variety of splitting 
schemes, which involve interactions of the re- 
latively large electronic magnetic moment with 
the field H, yield for H = 1000 oersted an anom- 
aly about twice the magnitude of the one we ob- 
serve in zero field. 

There remains the possibility that the energy 
levels of the nucleus are separated sufficiently 
to produce the anomaly. This splitting can be 
produced by interaction between the nuclear elec- 
tric quadrupole moment and the crystalline elec- 
tric field.* Because the nuclear magnetic mo- 
ment is much smaller than the electronic mo- 
ment, its interaction with the magnetic field is 
much smaller than in the electronic case, and 
leads to a relatively weak dependence of the spe- 
cific heat on H. The nuclei of the two stable Re 
isotopes have a spin of 5/2 and an average quad- 
rupole moment Q = 2.7 10-74 cm?. In the pre- 
sence of an axially symmetric field gradient, 
whose maximum value we denote by eq, six de- 
generate energy levels of the Re nucleus will be 
split into three doubly degenerate pairs, * the 
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separation of the lower two being 3e7qQ/20 and 

of the upper two, 6e7gQ/20. This splitting 

scheme leads to a specific heat A,T-?=(14/9)Ra,*T~* 
in the high-temperature limit T>>a,, where a, 

has been written for 3e7¢Q/20. Comparing this 
result with the experimental value of A,, we ob- 

tain a, = 0.002 deg and the maximum field gradi- 

ent eq = 0.1410** esu. Since the ingot is poly- 
crystalline and contains impurities, the value of 

eq may differ from that of a perfect crystal. 





*Supported by a Signal Corps Contract. 

‘This ingot, furnished to us by Dr. M. P. Garfunkel, 
has a history similar to specimen 4c of J. K. Hulm 
and B. B. Goodman [Phys. Rev. 106, 659 (1957)], for 
which: T, = 1.81°K with a transition breadth of 0.05 
deg. The known impurities are Fe=0.0011%, Si<0.01%, 
and W~0.03%. We are grateful for the loan of this 
sample. 

*Bardeen, Cooper, and Schrieffer, Phys. Rev. 108, 
1175 (1957). 

‘See, for instance, M. H. Cohen and F. Reif, Solid 
State Physics, edited by F. Seitz and D. Turnbull 
(Academic Press, New York, 1957), Vol. 5, p. 321. 

‘T, P. Das and E. L. Hahn, Solid State Physics, 
edited by F. Seitz and D. Turnbull (Academic Press, 
New York, 1957), Supplement 1. 








NUCLEAR SPECIFIC HEAT OF 
GALLIUM AND ZINC 


G. Seidel, * 
Kamerlingh Onnes Laboratorium, 
University of Leiden, Leiden, 
The Netherlands, 


and 


P. H. Keesom 
Purdue University, Lafayette, Indiana 
(Received February 25, 1959) 


Recently the specific heats of gallium’ and 
zinc’-* have been measured in the superconduct- 
ing state below 1°K. However the data have all 
been analyzed without concern of the possible 
effect of an interaction of the nuclear quadrupole 
moment with the gradient of the crystalline field. 
Now that such an effect has been shown to be of 
considerable importance for the specific heat of 
thenium,* it becomes necessary to review the 
previous results. 

Knight et al.> measured the pure quadrupole 
resonance spectra for Ga®® and Ga™ in the metal- 
lic state. From their results it becomes possi- 
ble to calculate the quadrupole coupling constant 


for both isotopes and also the contribution of the 
nuclei to the specific heat of naturally occurring 
gallium (60.2% Ga®*, 39.8% Ga”). This contribu- 
tion is given by C =4.6/T* ergs/mole deg when T 
>>0.001°K. It is only 5% of the specific heat of 
the conduction electrons at the lowest tempera- 
tures measured by us. Therefore we did not ob- 
serve a departure of the electronic specific heat 
in the superconducting state from the expected 
exponential decrease with temperature: 


Cog =¥T ea exp(-6T,/T). (1) 


Recalculating our data, after subtracting this 
nuclear term, the constants for Ga in Eq. (1) 
change: a from 7.0 to 7.5 and b from 1.35 to 
1.38. The fit of the experimental data and Eq. 

(1) is extended to somewhat higher temperatures: 
T,/ zh. 

Zinc has a hexagonal crystal, and hence an ap- 
preciable field gradient can exist at the nuclei. 
However the single isotope which is magnetic, 
Zn®” (spin 5/2), has a natural abundance of only 
4%. Phillips? noted deviations of the specific 
heat from the exponential temperature dependence, 
Eq. (1), outside the experimental error. It is 
reasonable to assume that these deviations are 
nuclear in origin and that they have a tempera- 
ture dependence proportional to 1/7°. From his 
graph it follows that the additional term is then 
given by C~1.5/T7* erg/mole deg for naturally 
occurring zinc. Hence for Zn®", C~37/T* 
erg/mole deg and the nuclear quadrupole coupling 
constant for Zn®™ is approximately 70 Mc/sec. 

If the additional term is subtracted from the spe- 
cific heat data of Phillips, his constants of Eq. 
(1) change: a from 5.8 to 6.2 and b from 1.22 to 
1.24 (using his values of y and T,,)- We reported 
values a=6.4, and b=1.27, which are only 
slightly higher. They were obtained at higher 
temperatures where this additional term can be 
neglected. 

In the normal state, at the lowest temperatures 
Phillips notes a 5% deviation from the expected 
value of the heat capacity but allows the possibil- 
ity of errors in his thermometers. However, 
this value of the deviation is to be expected if 
the splitting of nuclear energy levels of the Zn” 
isotope contributes to the specific heat. It should 
be possible to substantiate our assumption by a 
quadrupole resonance experiment. 

It is difficult to obtain data from Zavaritskii’s* 
curves for zinc, but it appears that his results 
are not inconsistent with the above conclusions. 

Zavaritskii® uses the relation a=4 for the 
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constants in Eq. (1). This relation holds very 
well for all the elements’ for which these con- 
stants are known. 


We would like to use this opportunity to correct 


a statement in our previous article.’ Dr. D. E. 
Mapother pointed out that on cooling through the 
transition temperature in a magnetic field a 
superconducting shell does not expel flux lines. 
Therefore our measurements were performed in 
the earth’s magnetic field. 





* 
National Science Foundation Postdoctoral Fellow. 


1G. Seidel and P. H. Keesom, Phys. Rev. 112, 
1083 (1958). 

2N. E. Phillips, Phys. Rev. Lett. 1, 363 (1958). 

IN. V. Zavaritskii, J. Exptl. Theoret. Phys. 
(U.S.S.R.) 34, 1116 (1958)(translation: Soviet Phys. 
JETP 34, 773 (1958)]. 

4p, H. Keesom and C. A. Bryant, preceding Letter 
[Phys. Rev. Lett. 2,260(1959)]. 

5 Knight, Hewitt, and Pomerantz, Phys. Rev. 104, 
271 (1956). 





THREE-LEVEL MASERS AS HEAT ENGINES* 


H. E. D. Scovil and E. O. Schulz-DuBois 
Bell Telephone Laboratories, 
Murray Hill, New Jersey 
(Received January 16, 1959) 


The purpose of this note is to demonstrate that 
three-level masers'>? can be regarded as heat 
engines. The principal conceptual difference 
between these and conventional heat engines is 
that in the 3-level maser one is concerned with 
the discrete energy levels of a particle’s internal 
energy whereas in a conventional heat engine one 
is concerned with the continuous spectrum of 
energies associated with external motion of the 
working substance. In treating a 3-level maser 
as a prototype of heat engine, a particular ad- 
vantage is, in our opinion, the resulting con- 
ceptual simplicity. Especially, it is easily shown 
that the limiting efficiency of a 3-level maser is 
that of a Carnot engine. 

Consider the system shown in Fig. 1. A three- 
level system is assumed with all transitions 
allowed and with no appreciable relaxation pro- 
cesses. The usual 3-level maser terminology is 
introduced by correlating transition 1—3 with 
pump frequency v,, 1-+-2 with signal frequency 
Ve, and 2-—3 with idler frequency v;. Asa 
further convention, the length of each energy 
level line is drawn proportional to its population. 
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FIG. 1. Three-level system in thermal contact with 
two heat reservoirs. 


The levels 1 and 3 are in thermal contact, 
through a filter passing frequencies in the vicinity 
of v, and rejecting frequencies in the vicinity of 
vy; and V5) with a heat reservoir at temperature 
T,. The temperature is indicated in the figure 
by showing schematically the Boltzmann distribv- 
tion of this heat reservoir. Levels 2 and 3 are 
in thermal contact with a reservoir at a lower 
temperature T, through a filter which passes 
frequencies in the vicinity of v; but rejects those 
close to v, and Ve: 

Experimentally, the high-temperature reser- 
voir might be realized by a gas noise lamp and 
the filter by a wave guide cutting off the lower 
frequencies. For practical purposes, however, 
the single mode present in a wave guide does not 
provide good thermal contact. The assumed 
coupling situation to the low-temperature reser- 
voir, on the other hand, was closely approxi- 
mated by experimental conditions in some maser 
experiments.* There, the idler transition of the 
gadolinium three-level system was coupled, 
through spin-spin interaction at frequency v;, to 
a transition of the same frequency of cerium 
ions within the same crystal. Thus, through the 
resultant short spin-lattice relaxation time, 
good thermal contact to the lattice heat reser- 
voir at T, was established. 

In the system described, for each quantum 
hvy supplied by the hot reservoir, the energy 
hv; is passed to the cold reservoir. The smaller 
quantum hv, can be extracted at the signal tran- 
sition if maser action prevails, that is if n,/n,?!. 
Thus the efficiency of this idealized system in 
maser operation is 


"yu = V5/Vp- (1) 


From the Boltzmann factors involved, we find 


Nz _ Ne Ns - hv; [hv (2) 
mo = ember!) ex (ee) 
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After rearrangement, this becomes 


Ng hv, V; T, = T> : 
Sooo Ea, 





In this formula, one recognizes the maser effi- 


ciency yg and the efficiency of the Carnot cycle, 
".* (T, - T,)/T,. , (4) 


Using these, the condition for maser action is 


| Iv Ne (5) 


This may be regarded as another formulation of 
the second law of thermodynamics. Maser effi- 
ciency equals that of a Carnot engine if the signal 
transition is at the verge of inversion, n,-n,~-+0 
or T sig ~e, 

As any heat engine, this system should be re- 
versible so that it acts as a refrigerator. This 
is indeed the case. Suppose a quantum hv, is 
applied to the signal transition. It causes an 
ion to go from state 1 to 2. The ion may further 
jump to state 3 if the energy hv; is supplied by 
the cold reservoir. The cycle is finally com- 
pleted when the ion returns from state 3 to state 
1 while the energy hv, is communicated to the 
hot reservoir. In this process, energy is ex- 
tracted from the idler transition, that is from 
the cold reservoir, so that it is refrigerated. 
The scheme outlined here, however, requires 
the signal transition to be absorptive. Other- 
wise the first step, application of hv,;, to the 
signal transition, would not have been possible. 
Thus the refrigeration scheme is possible if 
n,/n,<1. Again, the limiting efficiency of the 


refrigerator is that of a Carnot engine and it is 


realized with n, -,~-0 or Tgjg- +, 

It seems probable at this time that generation 
of microwaves through thermal excitation by 
two temperatures will be possible experimen- 
tally. Such a scheme should be very attractive 
for high microwave signal frequencies. Refri- 
geration experiments, on the other hand, as 
applied to interacting nuclear and electronic spin 
systems have been suggested by the theoretical 
work of Overhauser.* A thermodynamical analy- 
sis of the Overhauser effect has been given by 
Brovetto and Cini® and Barker and Mencher.® 

Finally, we should like to point out that the 
possibility of treating masers as heat engines 
sets a fundamental distinction between these and 
parametric amplifiers. Three-level masers are 
capable of operating with noise-like excitation in 
all three transitions. Parametric amplification 








requires some phase coherence, that is, mono- 
chromatic excitation is necessary for at least 
one of the three frequencies involved. This 
statement excludes heat as the source of energy 
for parametric amplification. 





“This work is partially supported by the Signal Corps. 

iN, G. Basov and A, M. Prokhorov, J. Exptl. Theoret. 
Phys. U.S.S.R. 28, 249 (1955) (translation: Soviet 
Phys. JETP 1, 184 (1955)). 

2N. Bloembergen, Phys. Rev. 104, 324 (1956). 

*Scovil, Feher, and Seidel, Phys. Rev. 105, 762 
(1957). 

“A, W. Overhauser, Phys. Rev. 92, 411 (1953). 

°P, Brovetto and G. Cini, Nuovo cimento 11, 618 
(1954). 

SW. A. Barker and A. Mencher, Phys. Rev. 102, 
1023 (1956). 





GAMMA-RAY ACTIVATION OF CARBON* 


L. D. Cohen and W. E. Stephens 
University of Pennsylvania, 
Philadelphia, Pennsylvania 

(Received February 26, 1959) 


A direct measurement of photonuclear reaction 
cross sections can be accomplished by the use of 
the monochromatic capture gamma rays of the 
T(p, y)He* reaction. 

If the gamma-ray energy be varied by changing 
the energy of the captured particle, the photo- 
nuclear cross section can be determined as a 
function of photon energy. The capture of pro- 
tons in tritium produces a high-energy gamma 
ray of (19.82 + {£p)-Mev energy, where Ey is the 
proton energy. This gamma ray has a smooth 
yield curve with no resonances, but it is of 
rather low intensity. Nevertheless, by using a 
highly effective detector we have measured the 
direct photoactivation of carbon by the C’*(y,n)C" 
reaction in the photon energy range 20.2 to 21 
Mev. 

The tritium was absorbed in a thin layer of 
zirconium evaporated on a silver backing’ which 
was water cooled. The thickness of the target 
for 1-Mev protons was measured to be 70 kev by 
observing the T(p,) threshold and comparing 
with very-thin-target differential curves.” 

The carbon was reactor-grade graphite in the 
form of an annular ring surrounding the tritium 
target. The carbon subtended the angles from 
90° to 102° with respect to the proton direction 
which corresponds to a range of photon energies 
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of about 50 kev due to the Doppler effect. These 
effects produced a spread in gamma-ray ener- 
gies, as shown in the upper part of Fig. 1, of 
about 81 kev at 20.2 Mev, 67 kev at 20.5 Mev, 
and 76 kev at 21 Mev. 

After a 40-minute irradiation at a given photon 
energy, the carbon disk was inserted in a posi- 
tron detector and counted for 20 minutes. Each 
point in Fig. 1 was the result of about 10 such 
runs. The positron detector consisted of two 
sodium iodide crystals 5 in. in diameter by 2 in. 
thick placed 3 in. apart, each viewed by a 6364 
Dumont photomultiplier tube, the pulses ampli- 
fied, channeled between 380 and 600 kev, and 
put into coincidence. This arrangement, with an 
18-in. iron and lead shield around it, had a back- 
ground of 3 counts per minute and an efficiency 
for positrons of 10% at the center tapering off 
smoothly to zero at the edge. 

The 20-Mev gamma rays were monitored by a 
large sodium iodide crystal 3 in. in diameter by 
4 in. thick. Pulses greater than 16 Mev and less 
than 22 Mev were counted. The amplifier and 
discriminator were calibrated on every run by 
recording the pulse-height distribution. 

The results of our measurements are given on 
the graph of Fig. 1. It is noticeable that no sharp 
resonances of the sort reported from betatron 
photoneutron yield curves® are evident in the en- 
ergy region examined. This is consistent with 
the total absorption measurements made in this 
laboratory on carbon with these monochromatic 
gamma rays.‘ It is also consistent with self- 
absorption measurements’ in carbon. 
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FIG. 1. Cy,m) activation curve. Cross section 
in millibarns for the photoproduction of C" radio- 
activity as a function of photon energy in Mev. The 
dashed figures indicate energy resolution vs energy 
deduced from tritium-zirconium target thickness 
measurements and Doppler angle calculations. 
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We would conclude from these measurements 
that the levels excited by photon absorption at 
this excitation energy (presumably mostly 17 
levels) have a level width approximately equal to 
their spacing or a spacing comparable with our 
resolution, or they could have an integrated 
cross section less than 0.05 Mev-mb. 

The help of Professor John Winhold in the ini- 
tial stages of planning this experiment is grate- 
fully acknowledged. 
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NUCLEON-NUCLEON SPIN-ORBIT 
INTERACTION POTENTIAL* 


M. H. Hull, Jr., K. D. Pyatt, Jr., 
C. R. Fischer, and G. Breit 
Yale University, 

New Haven, Connecticut 
(Received February 2, 1959) 


A recent publication’ reports the effect of a 
change in the spin-orbit interaction used by Signell 
and Marshak? and compares the results with the 
phenomenologic fits of Stapp, Ypsilantis and 
Metropolis.’ Calculations with the shortened 
range of the spin-orbit potential have been avail- 
able at Yale for some time but have not been pub- 
lished in view of questionable agreement with ex- 
periment. Additional calculation confirms the 
view that the representation of data by the poten- 
tial is in many respects not altogether satisfactory 
and indicates that the new fits to 310-Mev data’ do 
not appear reasonable in the light of the phase- 
shift energy dependence between 150 and 310 Mev 
and of disagreement with experiment within sev- 
eral times the experimental error in several cas¢s 
between 142 and 345 Mev. 

In Fig. 1 are shown phase shifts for the short- 
ened range of spin-orbit interaction denoted as 
MS 1 and for the original Marshak -Signell case, 
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curves) and *H, state (solid curves) as a function of 
energy calculated from MS, MS 1, and X potentials. 
For comparison, the SYM values of the °F. 2 Phase shift 
(circles) and °H, phase shift (squares) are plotted at 
312 Mev for their three best fits, labeled SYM 1,2,3. 


referred to as MS. Graphs employing a longer 
range* (16%) of the central singlet-even potential 
with the spin-orbit potential used only in P states® 
are shown and labeled by X. The three best fits 
of SYM are referred to by 1,2,3 starting with the 
best. The combination of F, and H, is not in agree- 
ment with SYM 1,2 or 3 and the claimed improve- 
ment in going from MS to MS 1 regarding agree- 
ment with SYM 1 does not lead to natural asso- 
ciation of the SYM 1 point with the theoretical 
curves. For some other phase shifts the agree- 
ment with SYM 1 is much better, notably so for 
%, *Py.1, ‘Dz. For some of these the X poten- 
tial gives the better agreement. 

At 120 Mev, MS 1 and X bracket the curve of 
polarization P versus center-of-mass angle 6 for 





$=50° to 90°, MS 1 is the better at 40°, and X is 
the better from 20° to 30°. However, for the 
ross section o the MS 1 values are too high by 
several times the experimental error while X 
represents the data appreciably better. At 142 
Mev there are disagreements with data outside 
experimental limits with both versions of the 


potential for o, while Po is represented well by 
both. The depolarization D is apparently rep- 
resented better by MS 1 than by X. At 240 Mev 
and 80<6@<90°, MS 1 differs from experiment for 
o by many times the statistical error while X 
does not disagree with experiment as badly; at 
intermediate angles the agreement is better for 
both MS 1 and X, while for @< 20° they agree with 
each other but disagree with experiment. A more 
complete account of the comparisons at this and 
other energies is being prepared for publication. 
There are many instances of radical departure 
from experiment both below and above 150 Mev. 
An example of a better fit by means of X to-D at 
312 Mev is shown in Fig. 2. In view of discrep- 
ancies between X and measurements in other re- 
spects, no special claims for X are made. It is 
rather desired to emphasize that both potentials 
have many deficiencies. A suitable velocity de- 
pendence of the potential can produce agreement 
with experiment but, without restrictions based 
on theoretical plausibility, such a fit has little 
significance. 

A grant of IBM-704 machine time by the Inter- 
national Business Machine Corporation is grate- 
fully acknowledged. 
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search, U. S. Army. 
‘Signell, Zinn, and Marshak, Phys. Rev. Lett. 1, 


416 (1958). A careful reading of this reference leaves 
no doubt regarding the recognition by its authors of the 
existence of serious limitations on the potentials dis- 
cussed by them. 

2p. S. Signell and R. E. Marshak, Phys. Rev. 109, 
1229 (1958). 

3Stapp, Ypsilantis, and Metropolis, Phys. Rev. 105, 
302 (1957), referred to in text as SYM. 

‘Fischer, Pyatt, Hull, and Breit, Bull. Am. Phys. 
Soc. Ser. Il, 3, 183 (1958). 

5By introducing the orbital angular momentum oper- 
ator explicitly in the interaction Hamiltonian, one can 
arrange for the vanishing of all but the L =1 interaction 
effects. In a certain sense, therefore, the curves X 
refer to a “potential.” 





EVIDENCE FOR THE A/=3 RULE* 


Frank S. Crawford, Jr., Marcello Cresti,T 
Roger L. Douglass, Myron L. Good, 
George R. Kalbfleisch, M. Lynn Stevenson, 

and Harold K. Tichot 
Lawrence Radiation Laboratory, 
University of California, 
Berkeley, California 
(Received February 16, 1959) 


It has been clear for some time that the ex- 
perimental data relating to strange-particle de- 
cays are suggestively close to the values predic- 
ted by the A/=3 selection rule.'’? The largest 
discrepancy has been between the limits*>? 


0.28 < B<0.38, (1) 


predicted for the neutral branching ratio B 

= P(K,°=1° + 2°)/[P(K,°—1° + 7°) + P(K,°=-7* +77)], and 
the values B =0.14+0.06 found in the Columbia 
propane bubble-chamber experiment,‘ and B 

= 0.06 (one event) in the MIT iron-plate cloud- 
chamber experiment.® The prediction (1) is ob- 
tained, instead of the prediction B=1/3 of the 
“pure” AJ=3 rule, if one assumes enough A/ 

= 3/2 to account for the decay K*~2r. 

The particle-mixture theory of Gell-Mann and 
Pais predicts, for the fraction f of short-lived 
(K,°) decays, 

f =K,°/(all K®) =3. (2) 
This result is expected from CPT invariance 
alone, i.e., CP invariance is not needed.® If, as 
seems likely, the 27 modes constitute practically 


all the short-lived decays, * and K,°~27 is negli- 
gible,” Eqs. (1) and (2) can be combined to pre- 
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dict, for AJ=3, 
0.31 <Ry < 0.36, (3) 


where Ry =P(K,°=n*+n )/(allK®). 


The corresponding experimental results (pre- 
vious to the present experiment) are Ry = 0.42 
+ 0.05 (Columbia), * and Ry =0.46 + 0.08 (Michigan 
propane chamber).® These are only in fair agree- 
ment with the prediction (3), but when combined 
with Columbia’s B/2=0.07 give f=0.49+0.08, in 
good agreement with the particle- mixture pre- 
diction (2). 

At the time of the CERN conference, we pre- 
sented preliminary results, based on 450 decay 
events observed in our hydrogen bubble chamber. 
for Rx; and for the A branching ratio R A The 
Al =4 rule predicts® the value 


R, =P(A=p+77)/(all A) =2/3. (4) 


Our preliminary values were in good agreement 
with the predictions (4) and (3). We reported no 
results on the K,°~7°+7° mode at that time and 
thus could not check the predictions (1) and (2). 

We have now completed our analysis, and re- 
port on a total of 1091 events. Our data show no 
contradictions with the predictions of the AJ=} 
rule. On the contrary they are in remarkable 
agreement with them. 

We have observed 227 double events (D = 227) 
in which both the A and the K° decay within the 
prescribed “fiducial volume” inside the chamber 
via their charged modes A~p+m~ and K,°=1* +1. 
There are 594 events in which only the A is ob- 
served to decay via its charged mode (A =594), 
and 270 events in which only the K,° is observed 
to decay via its charged mode (K =270). For 
each of the 864 single V’s, the production and 
decay dynamics of the observed particle checks 
with the hypothesis of associated production via 
either 7~ +p—A+K° (A production) or 1” +p->" 
+K° (=° production). 

The branching ratios Rx and Ra, and the “true 
number of associated productions, n(true), are 
given by 


” 


Rx =[D/(A+D)] nx, 
Ry =[D/(K+D)|n,, 
and 
n(true) =[(A+D)(K+D)/D]n, 


where NK» Na» and 7, are escape corrections 
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which would be unity for an infinite chamber. 

The n’s are calculated by averaging over 
smoothed “true” distributions in the production 
position and angle. We require that for a A or 
K° to be “detectable” its charged decay must oc- 
cur beyond 0.3 cm from the production point and 
must lie within the fiducial volume. The fiducial 
yolume is defined by the requirement that all de- 
cay tracks be at least 2 cm in length. In the 
small fraction of single V’s in which the produc- 
tion and decay dynamics of A and K° overlap 

when curvature and angle measurements alone 
are used, we measure ionization on the positive 
decay fragment, to distinguish protons (A decay) 
from 7* (K,° decay). 

The calculated averages of the n’s vary by only 
afew percent over the entire incident-pion mo- 
mentum range, and differ by only a few percent 
between the A and £° production modes. There- 
fore we quote here only the “grand average” 
over all the “true” distributions, which yields 
the calculated values Nk = 1.249, NA = 1.398, and 
j)=1.097. In performing the averages we use 
our lifetime values Ty =0.94+0.05 10°" sec 
and 7) =2.72+0.16x10"'° sec. We also calculate 
the derivatives with respect to lifetimes, with 
the results AlnRy =+0.193Aln7, and AlnR, 
=+0.149Aln7,. The contributions of uncertain- 
ties in lifetimes to the uncertainties in Ry and 
Ry are negligible. 


The observed counts, D, A, and K must be re- 
placed by D/e,, A/e,, and K/e,, where €, and «€, 
are scanning efficiencies for finding double and 
single V’s. By rescanning 30% of the film we 
find for “detectable” decays, ¢€, =0.995+0.005, 
and €, = 0.976 + 0.008. 

The values of Ry and Ra obtained at the var- 
ious beam momenta, 0.95 and 1.03 Bev/c (below 
t’ threshold), and 1.09, 1.12, and 1.23 Bev/c 
fabove 5° threshold) agree within the errors. 
Similarly the momentum-averaged A -production 
results agree within one-third standard devia- 
tion with the °-production results. We there- 
fore present only the grand average results, for 
the fractions of charged decays, and total num- 
ber of associated productions, 


Ry = 0.339 + 0.020, (5) 
R , = 0.627 + 0.031, (6) 
n(true) = 2020+ 100. (7) 


We turn now to the decay K,°-7°+7°. We have 
seen three events consistent with this decay 








mode; in each case there is an associated 
charged A decay. In one case a 7° undergoes a 
“Dalitz decay” into e++e°++y. In this other two 
cases one of the 7°-y rays produces an electron 
pair in the liquid hydrogen. From our experi- 
mental K® momentum distribution we find the 
7°-y-ray spectrum by assuming isotropy in the 
K° and n° decay. Combining this with the known 
pair-production cross section per hydrogen 
atom, and including a probability for Dalitz-de- 
cay of 2.510 per K,°, we find a total detection 
efficiency per K,° of 3.5x107*. Our three events 
thus correspond to 3/3.5 x10~? = 86+ 50 decays. 
The number of accidental counts due to chance 
coincidences from unassociated electron pairs is 
estimated from the frequency of pairs, and the 
chance of fitting the decay dynamics. The result 
is that less than 0.2 accidental count is expected. 
(No correction was made.) There were 227 
decays K,°=1*+7~, associated with charged A 
decays. Therefore, independent of assumptions 
as to the value of f=K,°/(all K°), and independ- 
ent of escape corrections, we find for the frac- 
tion B of neutral K,° decays 


B =86/(86 + 227) =0.27+0.11, (8) 


which is consistent with prediction (1). We com- 
bine the charged and neutral results (5) and (8) 
to obtain 


f=0.47+ 0.080, (9) 


in good agreement with the particle-mixture- 
theory prediction (2). 

The disagreement between our results (5), (8), 
and (9) and those of other groups may be ex- 
pressed as follows. Both Eisler et al.* and we 
find the “expected” f=4, but disagree on the ad- 
mixtures of charged and neutral K,° decays. 
However, if we combine all experiments without 
imposing the constraint f= then our charged K,° 
decays dominate statistically, whereas in the 
neutral decays the Columbia and MIT results 
prevail (7+1-=8 events to our 3 events). The in- 
ternal consistency is furthermore not bad—a 
probability of 18% for the charged decays, and 
20% for the neutral. The resulting charged and 
neutral world averages are R,(U.C., Col.,* 
Mich.®) =0.354+0.018, and B(U.C., Col.,* MIT®) 
=0.141+0.037. (The error on B reflects only the 
counting statistics.) In the absence of other de- 
cay modes of the (short-lived) K,° these combine 
to give f=0.424+0.026, in disagreement with 
prediction (2). If, instead, f=5 is assumed, 
then 18 + 6% of the K,° decays would have to occur 
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by an as yet unknown mode. We believe that it 
is instead more reasonable not to combir the 
experiments, and to attribute the disagreement 
to statistical bad luck or systematic error. The 
¥ probability that our charged fraction is con- 
sistent with the Columbia-MIT neutral fraction 
is then about 107°. 

We next consider the decay A-n+7°. We have 
found two Dalitz decays and one y conversion 
corresponding, presumably, to this mode. In 
each case there was an associated charged K,° 
decay. The three events correspond to 171 +100 
neutral A decays. Combining these with the 227 
double V’s, we find R, =227/(227 +171) =0.57 
+0.14. Since we have good reason for believing 
that the associated-production hypothesis is 
valid and that there are no prominent A decay 
modes other than the two considered here, this 
last result can be combined with result (6) for 
the charged decays. We thus obtain the weighted 
average 


Ry= 0.624 + 0.030. (10) 


Result (10) is in excellent agreement with the 
Al=4 prediction R =2/3, and with other deter- 
minations by groups at Columbia, *~* Michigan, ® 
and MIT,°*~® as well as with results of the Berke- 
ley K -capture experiment.® A least-squares 
weighted average of these results gives R AtU.C., 
Col., Mich., MIT) =0.637 + 0.020. 

Finally, Dalitz' and Pais and Treiman’® point 
out that if the AJ=4 rule is valid one expects 
w(K,°-1* + 7° +7°) = #w(K*=37). Furthermore, 
the decay rate w(K,°~1*++7~ +7°) should be ex- 
ceedingly small compared with the correspond- 
ing K,° rate,'° so that we may attribute any ob- 
served 1*+7~+7° decay to the K,°. 

We have seen one such decay. The event was 
associated with a charged A decay. However, 
we would easily detect this decay mode if it oc- 
curred as a single V. Corresponding to the total 
number of associated productions, n(true) =2020, 
there should be 1010 K,°’s. Combining the 
known K* branching ratios" and lifetimes’? and 
the AJ=4 rule, one obtains the prediction 
w(K,°~1* +7" + 7°) =2.410° sec™'. The average 
potential time for K*’s was 3.5x107'° sec. We 
therefore expect to find 2.4x10°®x3.5x10 % 
x1010 = 0.85+0.10 decays K,°-1*+7>+7°. Our 
single event is thus consistent with the A/J=3 
rule.*® 

In summary, we find that our 1091 associated- 
production and decay events are in remarkably 
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good agreement with the predictions of the a/- 
rule. 

This makes the total experimental evidence 
for AJ =3 fairly impressive: (a) The decay 
= -n+mn~ occurs, if at all, much less frequently 
than = ~A+7~. This can be understood if the 
= has J=4 and the AJ=4 rule holds. (b) The 
strong inhibition of K*~27 relative to K,°~27 
follows from the AJ=4 rule, if the K has zero 
spin and /J=3. (c) The admixture of A/J=# re- 
quired to admit the observed K*~2z rate is in 


4 
2 





good agreement with our K,°~27 branching ratios, 


(d) The A branching ratio agrees with (but does 
not require®) AJ=3. (e) The branching ratio’ 
P(K*=1* 407° +2°)/P(Kt=n* + 27°) agrees with (bu 
does not require’) AJ=3. (f) The results of 
Cool et al.’® on =* decay asymmetry are most 
easily explained by (but do not require) A/=}. 
(g) Our one K,°=7* +7 +7° agrees (as well as 
one event can) with AJ=3.'° 

We wish to thank Luis W. Alvarez for his in- 
terest and support, and Don Gow and Hugh 
Bradner for their assistance. 
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Energy Commission. 
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CHARGE INDEPENDENCE IN THE REACTIONS 
p+d=n°+He® AND p+d—1*++H® AT 450 Mev* 


A. V. Crewe, E. Garwin, B. Ledley, 
E. Lillethun, R. March, and S. Marcowitz 
Enrico Fermi Institute for Nuclear Studies, 
University of Chicago, 

Chicago, Illinois 
(Received February 19, 1959) 


According to the principle of charge independ- 
ence, the ratios of the cross sections for the 
two reactions p+d—n++H®* and p+d—n°+ He® 
should be exactly a factor of two. This statement 
should be true both for the total cross section 
and the differential cross section. A previous 
attempt to establish this ratio has been made’ in 
which measurements of E and dE/dx were used 
to identify the He*. Owing to the large back- 
ground of protons and deuterons, this separation 
was difficult to achieve. In order to measure 
the cross section for the second reaction, both 
the »*+ and the H® in coincidence were detected. 

In the present experiment the differential cross 
sections for these two reactions were measured 
by counting only the heavy particles. First of 
alla momentum analysis was performed in order 
to substantially reduce the background. This 
momentum analysis was performed with the aid 
fa magnet which has been previously described.” 
The deflected particles were counted with a quad- 
tuple coincidence arrangement consisting of three 
counters in coincidence and the last one in anti- 
coincidence. An absorber was placed between 
the last of the three counters and the anticoinci- 














DIFFERENTIAL CROSS SECTION 
ARBITRARY UNITS 






dence counter. Upon the receipt of a signal show- 
ing that a particle had stopped in the absorber 

the pulse height in the third counter was meas- 
ured with a 50-channel pulse-height analyzer. 
This measurement of pulse height served to dis- 
tinguish the H* and He® from the background of 
protons and deuterons. 

The external proton beam from the University 
of Chicago’s synchrocyclotron was focused into 
a spot about 1/2 in. in diameter on a polymerized 
CD, target. The thickness of this target was 
chosen to be approximately 0.8 mm. This was 
necessary because the range of the lowest mo- 
mentum He® particle is only 4 mm of plastic 
scintillator. The secondary particles from the 
target were then deflected by the magnet into the 
counter telescope. The whole system of proton 
beam pipe, target, and magnet was maintained 
under vacuum at a pressure of a few microns. 
The experiment consists of taking CD,-C differ- 
ences at various magnet current settings and 
measuring the pulse-height spectrum at each 
setting. The magnet can approach within 8 de- 
grees of the forward direction without the proton 
beam striking the magnet, and the maximum 
angle of emission of either of the heavy particles 
is 15 degrees. We have measured the production 
of these heavy particles from 8 to 13 degrees in 
the laboratory system for the low-momentum 
portion of the spectrum. Owing to the small 
range of the He®* particles, it was necessary to 
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FIG. 1. The differential cross sections for the reac- 


tions p+ d-—-1* + H® and p+d--n°+ He® with 450-Mev 
protons. The measured cross sections for the second 
reaction have been multiplied by a factor of two in order 
to facilitate a direct comparison. The cross sections 
are given in arbitrary units which are, however, ap- 
proximately microbarns/steradian. 
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use only three counters in this part of the ex- 
periment. However, measurements on the H* 
were repeated in order to show that the efficiency 
of detection was identical in the two cases. 

The differential cross section of these two pro- 
cesses at the angles measured is shown in ar- 
bitrary units in Fig. 1. The differential cross 
section for the production of He* has been multi- 
plied by a factor of 2 on this graph. A weighted 
mean for the branching ratio ogs/oy,3 has been 
computed from these data to be 1.91+0.25. Ac- 
cordingly we conclude that there is no evidence 























ee, 


to show that the principle of charge independenc, 
does not hold in these reactions. 

The authors would like to acknowledge the cop. 
siderable help given by Charles Rey in preparing 
for and carrying out these measurements. 
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In this section are printed the abstracts of Articles that 
have been forwarded to The American Institute of Physics 
for publication in THE PHYSICAL REVIEW. In quoting 
information obtained from this section before the appear- 
ance of the corresponding Article, reference should be 
made to “Physical Review (to be published)” rather than 
to this Journal. 





STABLE ORBITS OF CHARGED PARTICLES IN 
AN OSCILLATING ELECTROMAGNETIC FIELD. 
Erich S. Weibel, Physical Research Laboratory, 
Space Technology Laboratories, Incorporated, 
Division of The Ramo-Wooldridge Corporation, 
Los Angeles, California (Received November 3, 
1958). 


Orbits in a field of the cylindrical wave guide 
driven in the TE,, mode have been studied. The 
radial components of the Lorentz force acting on 
the particle is never positive regardless of its 
charge; thus, it is attracted toward the axis of 
the wave guide around which it moves in a com- 
plicated path. Sufficient conditions are given 
under which the orbits are stable, that is, remain 
bounded for all times. 


NATURE AND ROLE OF IONIZING POTENTIAL 
SPACE WAVES IN GLOW-TO-ARC TRANSI- 
TIONS. Russell George Westberg, The Depart- 
ment of Physics, The University of California, 
Berkeley, California (Received July 7, 1958). 


The transition from a glow discharge to a 
transient arc, induced in Al-cathode tubes 145 
cm long by shorting of a series resistor, has 
been studied for a variety of gases by observing 
currents at the electrodes, the potentials at 
various points in the tube using probes, and the 
movement of luminous pulses using photomulti- 
plier tubes and fast oscilloscopes. With a dc 
power supply with up to 2600 volts applied and a 
large capacity parallel to the discharge tube, 
aid with due regard to proper impedance match- 
ing, relatively pure gases N,, A, H,, 0,, and 
mixtures of N, and 0, were studied from 50 mi- 
trons of Hg up to several hundred microns. In 
udition, shorter tubes with Cu, Ni, Pt, Hg, and 
Weathodes and filled with He or N,, at interme- 
diate pressures up to 7.2 mm of Hg, were brief- 
ly studied. The decrease of the series resistance 
increased the very low current in the normal 
slow mode to higher values up to 50 ma usually 








in the abnormal mode. With slightly oxidized 
cathodes, there followed, after times varying 
from 10°° second to many seconds, catastroph- 
ic breakdowns to a transient power arc of the 
order of 10-20 or more amperes. The transi- 
tion was initiated by a burst of electrons from 
the cathode yielding a cathode current increase 
to several amperes in 2-10x10~* second. This 
burst sends a luminous pip or pulse to the anode 
at velocities varying from 5x10® to 8x10° cm/sec, 
depending on pressure, voltage, and gas. On 
arrival at the anode, a luminous return arc plas- 
ma front moves towards the cathode with speeds 
from 10° to 10° or more cm/sec. When it ar- 
rives at the cathode, the conductivity is suffi- 
ciently homogenized along the column so that the 
current and luminosity increase along the whole 
column until the drain on the capacity lowers the 
supply voltage. The luminous pip and, to a less 
clearly defined degree, the return plasma front 
are accompanied by potential space waves which 
can be established by cross-plots yielding lumi- 
nosity and potential distributions across the tube 
at various times. The transition appears to be 
initiated by a sudden breakdown of the oxide 
layer on the cathode as it thins by ion bombard- 
ment, causing a small fraction of the oxide mol- 
ecules to be ionized within millimicroseconds, 
thus yielding the initial current burst. These 
electrons, accelerated by nearly the full poten- 
tial drop across the dark space, move as ion- 
izing potential fronts, accompanied by luminos- 
ity, with velocities in most cases dictated by the 
cathode fall potential. In O,, when pressures 
permit effective photoionization within the tube 
length, the velocities of the front are increased 
up to 8x10° cm/sec. 


AC-DC ELECTROLUMINESCENCE. W. A. 
Thornton, Research Department, Westinghouse 
Electric Corporation, Bloomfield, New Jersey 
(Received April 24, 1958; revised manuscript 
received August 27, 1958). 


The enhancement of electroluminescence, by 
the superposition of ac and dc voltages, has 
been observed in certain electroluminescent 
zinc sulfide powder phosphors. The addition of 
a dc voltage to an ac voltage exciting visible 
electroluminescence increases the emission as 
much as 250 times, under conditions where the 
de electroluminescence alone is about equal to 
the initial ac electroluminescence. Wave-form, 
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transient, and spectral measurements show 
that the increased emission is recombination, 
at dc-ionized luminescence centers, forced by 
the ac field. 


MAGNETIC FIELD DEPENDENCE OF THE 
SURFACE IMPEDANCE OF SUPERCONDUCT- 
ING TIN. M. Spiewak, Department of Physics 
and Institute for the Study of Metals, University 
of Chicago, Chicago, Illinois (Received October 
27, 1958). 


The surface impedance at 1 kMc/sec of super- 
conducting Sn has been investigated as a func- 
tion of temperature, 7, and of static magnetic 
field, H, both longitudinal and transverse to the 
rf current. The following unexpected phenomena 
have been observed. The rf resistance, R, and 
the reactance, X, may each decrease with in- 
creasing H. A decreasing X and an increasing R 
may coexist. A variation of either the crystal- 
line or static field orientation with respect to 
the rf magnetic field can change the sign of 
[R(H)-R(0)] and of [x(H)-x(0)]. For T<2.0°K, 
both R and X are monotonically increasing func- 
tions of H. In this range of 7, the approximation 
of a quadratic dependence for R on longitudinal 
and transverse H becomes increasingly unsatis- 
factory with decreasing T. The quadratic coeffi- 
cient of the variation of X with longitudinal H 
also decreases with decreasing 7, but the cor- 
responding coefficient in transverse H increases, 
1.2°K<T<2.0°K. The relation between these un- 
expected behaviors and previous experimental 
and theoretical studies is discussed. 


EXCITATION SPECTRA AND TEMPERATURE 
DEPENDENCE OF THE LUMINESCENCE OF ZnS 
SINGLE CRYSTALS. A. Halperin and H. Arbell, 
Department of Physics, The Hebrew University, 
Jerusalem, Israel (Received September 3, 1958). 


The luminescence of ZnS:Cl and ZnS:Cu:Cl 
crystals was measured for the temperature re- 
gion 80-500°K and for different wavelengths of 
exciting light. The behavior of the luminescence 
versus temperature curves differed from similar 
curves for powders reported in the literature. 

Excitation spectra were also obtained for the 
green and blue luminescence bands. These were 
found to consist of a complexity of fine bands. 

Each of the bands was found not to shift when vary- 
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ing the temperature of the crystal. The relative 
intensity of the bands near the absorption edge oj 
the crystal, however, changed with temperature 
in such a way as to cause an apparent shift in the 
maximum of the excitation spectrum in this regi, 


“THIN- FILM” EXPERIMENT WITH BULK SU- 
PERCONDUCTORS. Hans Meissner and Richard 
Zdanis, The Johns Hopkins University, Balti- 
more, Maryland (Received October 9, 1958), 


Hollow cylindrical samples of indium have been 
prepared by an extrusion process and equipped 
with both a center wire and a concentric tube. 
The critical current has been determined for the 
case in which the current passes through the 
sample and returns through the concentric tube. 
It has been found that this critical current is 
identical with that of solid samples. The same 
determination has been made for the case in 
which the current returns through the center 
wire. In this case the critical current is about 
80 to 90% of that for solid samples. A field was 
produced by the center wire, the respective 
current returning through the concentric tube. 









The resistance of the sample was measured with 
a small measuring current, also returned 
through the concentric tube, as a function of the 
field produced by the center wire. The critical 
field thus determined depends on the value of 
the measuring current. These findings can be 
interpreted in the following way: If the magnetic 
field of the current through the center wire is 
dominating, the current in the sample will flow 
in a very thin layer at the outer surface of the 
sample. This constitutes a “thin-film experi- 
ment” which does not require a thin metal film. 
The resistance measurements have been supple- 
mented by measurements of the circular flux on 
a sample of larger size. 


EFFECT OF PRESSURE ON THE INFRARED 
ABSORPTION OF SEMICONDUCTORS. L. J. 





Neuringer,* Department of Physics, University 
of Pennsylvania, Philadelphia, Pennsylvania 
(Received May 26, 1958; revised manuscript re- 
ceived December 1, 1958). 


Measurements have been made of the effect of 
hydrostatic pressure upon the intrinsic infrared 
absorption of germanium, silicon, and tellurium 
in the pressure range 1-2000 atmospheres. In 
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germanium the variation with pressure of the 
lowest lying conduction band minimum, (111), 

js found to be (7.3+1.5)x10~* ev/atmos. It is 
suggested that the variation with pressure of the 
next higher conduction band edge, (000), hasa 
value of (11+3)x10~* ev/atmos. The shift of the 
infrared absorption edge with pressure in silicon 
is small and is toward smaller energy, amount- 
ing to about -2 x 10~° ev/atmos between 1 and 

9020 atmospheres. The sign and magnitude of 

the pressure coefficient in this material are in 
agreement with the results of Warschauer, Paul, 
and Brooks. For tellurium both Eg, the energy 
gap corresponding to light polarized perpendic - 
war to the c axis of the crystal, and E,,, the 
energy gap for polarization parallel to the c 

axis, decrease as the pressure is increased. It 
appears that E,, decreases more rapidly with 
pressure than EF, The mean pressure coeffi- 
cient of Eg, is ay. 2+0.4)x10~ ev/atmos. For 
Eg, the ane pressure coefficient is -(1.8+0.3) 
x10*% ev/atmos. The measured pressure co- 
efficients are used to calculate the thermal dila- 
tion term in the equation for the change of the 
energy gap with temperature for each material. 
The electron-lattice interaction term appearing 
in this equation is then deduced. In silicon these 
two terms are of opposite sign with the electron- 
lattice term dominant. In germanium the elec- 
tron-lattice interaction effect accounts for 75 
percent of the effect of temperature on the ener- 
gy gap. For tellurium the two effects are almost 
equal in magnitude but of opposite sign. 


"Now at Research Division, Raytheon Manufacturing 
Company, Waltham, Massachusetts. 


THEORY OF ULTRASONIC CYCLOTRON RESO- 
NANCE IN METALS AT LOW TEMPERATURES. 
1. Kjeldaas, Jr., Westinghouse Research Labo- 
ratories, Pittsburgh, Pennsylvania (Received 
August 25, 1958). 


The extent to which measurement of the propa- 
gation characteristics of circularly polarized 
acoustic shear waves as a function of a magnetic 
field applied parallel to the direction of propa- 
gation can reveal information about the band 
structure of metals is investigated theoretically. 
The interaction between the electrons and sound 
waves is treated classically and is viewed as 
taking place via an internal electric field and by 


scattering modified by the sound wave. The treat- 





ment is based on the combined solution of Boltz- 
mann’s and Maxwell’s equations for a free- 
electron model. A simple physical picture in- 
volving cyclotron resonance is developed by 
means of which the treatment is extended to 
more general band models. The most striking 
results are obtained in case the mean free path 
is larger than the wavelength. Under these con- 
ditions, it is found that subject to certain re- 
strictions, (1) absorption edges exist, the 
measurement of which can be related to the cur- 
vature of the Fermi surface; (2) the shape of the 
Fermi surface may be determinable from the 
shape of attenuation and dispersion curves within 
the absorption region; and (3) the cyclotron fre- 
quency may be measured. The results (2) and 
(3) may be dependent upon the simplifying as- 
sumption regarding the electron-sound wave 
interaction. 


ROLE OF DOUBLE -QUANTUM TRANSITIONS 
IN MASERS. Shaul Yatsiv, Stanford University, 
Stanford, California (Received July 7, 1958). 


A generalized treatment of a three-level maser 
as an example of double-irradiation experiments 
is based on a previous analysis of multiple-quan- 
tum transitions in nuclear magnetic resonance. 
The emphasis in the present treatment, as dis- 
tinct from other expositions of the same problem, 
is on the role of the double-quantum transitions. 
In particular, conditions are found in which the 
maser operation is mainly governed by the dou- 
ble-quantum process and does not require a true 
“pumping” stage. Such a case is realizable in 
practice but its technical applicability is ques- 
tionable. A study of the behavior of various tran- 
sitions in this example is instructive for under - 
standing the general role of the double-quantum 
transitions in double-resonance experiments. 


NUCLEAR STRUCTURE EFFECTS IN INTERNAL 
CONVERSION COEFFICIENTS BY CONFIGURA- 
TION MIXING. L. S. Kisslinger, Western Re- 
serve University, Cleveland, Ohio, and Oak Ridge 
National Laboratory, Oak Ridge, Tennessee (Re- 
ceived June 2, 1958). 


The nuclear matrix elements which are needed 
to determine the internal conversion coefficients 
when a finite nucleus is employed are derived for 
nuclei for which the shell-model wave functions 






273 








VOLUME 2, NUMBER 6 


PHYSICAL REVIEW 








LETTERS MARCH 15, 1959 





are a good zero-order approximation for low- 
energy processes. Using configuration mixing, 
general expressions are derived for these matrix 
elements. It is shown that the nuclear structure 


alteration can be ten to twenty percent or more 
for / -forbidden transitions. Numerical results 
are given for the M1 and E2 279-kev transitions 
in TF. 


SEARCH FOR POSITRON EMISSION IN K*. 

D. R. Tilley and Leon Madansky, Johns Hopkins 
University, Baltimore, Maryland (Received 
May 19, 1958). 


A search for positron emission in K* has been 
made with the use of a triple coincidence tech- 
nique. An upper limit of (3.6+1.8) x10™* posi- 
tron per second per gram of natural potassium 
has been set. A corresponding upper limit of 
0.57+0.29 has been computed for the ratio of the 
squares of the matrix elements, M jr.” , for 

the K* g* and 8 transitions. 


ANALYSIS OF THE REACTION C*“(He’, a)C’?. 
George E. Owen, L. Madansky, and Steve 
Edwards, Jr., Department of Physics, The Johns 
Hopkins University, Baltimore, Maryland (Re- 
ceived July 22, 1958). 


The angular distributions of the alpha particles 
from the reaction C'"(He*, a)C’” at 4.5 and 2.0 
Mev have been analyzed from the point of view of 
nuclear stripping. By symmetrizing the wave 
function of the complete Hamiltonian with re- 
spect to the exchange of alpha particles, the 
heavy-particle stripping amplitude is formally 
included. The experimental distributions change 
decidedly between 2.0 and 4.5 Mev, and the 
calculated curves employing a combination of 
pick-up and a-stripping from C** agree quite 
well with the experimental distributions. Inter- 
ference between the two channels plays a major 
role in the final expressions. 


QUANTUM-MECHANICAL THREE-BODY PROB- 
LEM. Leonard Eyges, Lincoln Laboratory, 
Massachusetts Institute of Technology, Lexing- 
ton, Massachusetts (Received November 3, 1958). 


We treat the quantum-mechanical problem of 
three spinless particles, with the boundary con- 
dition that the logarithmic derivative of the wave 
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function be a prescribed constant at each of the 
three boundaries Ir, -r,|=a, ir, -T,! =a,IT, “1, 
=a. This boundary condition is discussed; it is 
roughly equivalent to an interparticle potentia] 
which consists of a hard core plus a strong shor. 
range attractive part. The eigenfunctions and 
eigenvalues of the system are given by the solu- 
tions of an infinite set of coupled homogeneous 
integral equations. The equations involve partia| 
wave expansions in the interparticle distances 
and can often be truncated with good approxima- 
tion by taking only a finite number of partial 
waves. We discuss the solution of these equa- 
tions for the ground state of the system, taking 
relative S-waves only, for which case the infinite 
set of equations reduces to a single integral 
equation in one variable. 


ALPHA-DECAY BARRIER PENETRABILITIES 


WITH AN EXPONENTIAL NUCLEAR POTENTIAL; 


EVEN-EVEN NUCLEI. John O. Rasmussen, Ra- 
diation Laboratory and Department of Chemistry, 
University of California, Berkeley, California 
(Received October 27, 1958). 


The real potential derived by optical -model 
analysis of data on elastic scattering of alpha 
particles is used for calculation of barrier pen- 
etrabilities for all known alpha decay groups of 
even-even nuclei. The barrier penetration factors 
were calculated by numerical integration in the 
WKB approximation taking into account centrif- 
ugal barrier effects but ignoring noncentral in- 
teractions. Using these penetration factors and 
the experimental alpha half-lives, the reduced 
level widths 5? are calculated. Ratios of 6? values 
for ground and excited-state alpha groups are 
tabulated as a set of reduced hindrance factors. 


RATIO OF ASYMMETRIC TO SYMMETRIC 
FISSION OF U?* AS A FUNCTION OF NEUTRON 
ENERGY. R. B. Regier, W. H. Burgus, and R. L. 
Tromp, Phillips Petroleum Company, Atomic 
Energy Division, Idaho Falls, Idaho (Received 
October 15, 1958). 


Radiochemical measurements of the relative 
yields of the fission products Mo”, Ag'!!, Ag”, 
and Cd'!5 from the low-energy neutron-induced 









~~ 








fission of U2** have been made at various neutron | 


energies. The energies chosen were thermal, 
1.8, 2.3, and 4.7 ev, the latter three correspond- 
ing to the peaks of previously reported resonances 
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It was found that the ratio of asymmetric to sym- 
metric fission is larger at the 1.8- and 2.3-ev 
resonances than at thermal energies. At the 4.7- 
ev resonance, however, this ratio is the same 

as at thermal energies, to within experimental 
uncertainties. In addition, it was found that the 
ratio for epi-cadmium neutrons differed from 

that for thermal neutrons. The results are con- 
sistent with Wheeler’s prediction that the ratio 

of asymmetric to symmetric fission should de- 
pend upon the spin state of the fissioning nucleus. 


DECAY OF As™. D. J. Horen, W. E. Meyerhof, 
J. J. Kraushaar,* D. O. Wells, E. Brun,f and 

J. E. Neighbor,} Stanford University, Stanford, 
California (Received October 13, 1958). 


The decay of As” (17.5 day) has been reinvesti- 
gated with single and coincidence scintillation 
spectrometers. The singles spectrum revealed 
previously unknown gamma rays of 1.19+0.01, 
1.60+0.04, and 2.22+0.02 Mev. X-ray-gamma 
coincidence measurements indicate the 1.19-Mev 
gamma ray is predominantly associated with the 
decay to Ge™*. Evidence is found for 0.6—0.6- Mev 
gamma-gamma coincidences in Ge“. Based 
on these data, a level scheme for Ge” is pre- 
sented with excited states at 0.596, 1.19, and 
2.22 Mev. No new information concerning the 
decay to Se™* has been gained. 


"Present address: Department of Physics, University 
of Colorado, Boulder, Colorado. 

‘Present address: Physics Department, University 
of Zurich, Zurich, Switzerland. 

Present address: Massachusetts Institute of Tech- 
nology, Cambridge, Massachusetts. 


LOW-INTENSITY CONVERSION LINES FROM 
Sm‘** AND Sn¥5, W. T. Achor, Vanderbilt Uni- 
versity, Nashville, Tennessee, and Western Re- 
serve University, Cleveland, Ohio, and W. E. 
Phillips, * a. ¥. Hopkins, f and S. K. Haynes, t 
Vanderbilt University, Nashville, Tennessee 
(Received November 13, 1958). 


The electron spectrum of Sm!*! has been meas- 
ured with a lens spectrometer and its gamma- 
Tay spectrum with a scintillation spectrometer. 
The L-conversion coefficient of the 21.7+0.3 kev 
de-excitation of the first excited state of Eu'* 
is 2024, with a7 /ays,y =2.2+0.4 and ay/ay 
=2.2+0.2, indicating M1 multipolarity. The low 


intensity of the conversion electrons relative to 
the continuous beta spectrum of end point 75.9 
+0.6 kev implies the existence of a weak beta 
transition whose end point is 54.2+0.7 kev. The 
beta branching ratio, 54.2-kev beta transitions / 
75.9-kev beta transitions, is 1.7/100. 

Log,,ft values for the two beta decays are 
7.6+0.2 for the 75.9-kev decay and 8.8+0.2 for 
the 54.2-kev decay. A partial analysis of the de- 
cay in terms of the Nilsson-Gottfried scheme has 
been made, but no definitive assignment for the 
21.7-kev state appears to be possible. Measure- 
ments of the electron spectrum of Sn’* with the 
lens spectrometer reveal the existence of a weak 
conversion line corresponding to a transition of 
253 kev. The intensity of this conversion line is 
(2.8+1.0) x10 of the intensity of the conversion 
line from the 392-kev transition. This relative 
conversion electron intensity indicates that the 
253 -kev transition leads to the 392-kev state of 
In™*, implying the existence of an excited state 
at 645 kev. The existing uncertainty in the rela- 
tive gamma-ray intensities for the two transitions 
makes it impossible to assign spin and parity to 
the 645-kev state; the decay scheme presented 
is therefore incomplete. 


*Now at Texas Technological College, Lubbock, Texas. 
Now at Davidson College, Davidson, North Carolina. 
TNow at Michigan State University, East Lansing, 
Michigan. 


H°(a,y )Li’ AND He*(a, y )Be’7 REACTIONS. H. D. 
Holmgren and R. L. Johnston, Nucleonics Division, 
United States Naval Research Laboratory, Wash- 
ington, D. C. (Received October 14, 1958). 


The cross sections for the H*(a, y)Li’ and 
He*(a, y)Be’ reactions have been measured at 
bombarding energies of 0.48, 0.72, 0.94, 1.13, 
and 1.32 Mev. The cross sections at the lower 
bombarding energies may be fitted within the 
experimental uncertainties by the expressions 


o =0.12(1-0.51 Eg)Eqg™ exp(-3.96Eg~2) ub 
for the H® reaction and 
o =2.8(1-0.55Eq)Ea™ exp(-7.92Eq~2) ub 


for the He® reaction. 


HEAVY -MESON—HYPERON RELATIVE PARITY 
AND THE SIGN OF THE COULOMB INTERFER- 
ENCE IN K -p SCATTERING. E. Galzenati and 
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B. Vitale, Istituto di Fisica Teorica dell’ Univer- 
sita, Napoli, Italy, and Scuola di Perfeziona- 
mento in Fisica Teorica e Nucleare del C.N.R.N., 
Napoli, Italy (Received October 31, 1958). 


Two subtracted dispersion relations are ap- 
plied at the same time to the study of the depend- 
ence of the real part of the K -p scattering am- 
plitude at zero energy, D_(K), on the values of 
the coupling constants g,,”. The results, when 
compared with the existing experimental infor- 
mation, strongly indicate that the K meson is 
ps-ps (namely pseudoscalar with respect to both 
the A and the = hyperon) and that the sign of 
D_(w) at low energy is positive, giving therefore 
constructive Coulomb interference in the K -p 
elastic scattering. 


EFFECTIVE MEAN FREE PATH OF ANTINUC- 
LEONS IN NUCLEAR MATTER. Jose R. Fulco, 
Radiation Laboratory, University of California, 
Berkeley, California ( Received November 6, 
1958). 


The effective mean free path for antinucleons in 
nuclear matter is calculated at intermediate en- 
ergies from the Fermi gas model, in terms of 
the nucleon-antinucleon total and differential scat- 
tering cross sections. The results are used to 
obtain the imaginary part of the optical- model 
potential in the formalism of Riesenfeld and 
Watson. 


DECAY OF ,H*® AND THE SPIN DEPENDENCE 
OF THE A-NUCLEON INTERACTION. M. Leon, 
Laboratory of Nuclear Studies, Cornell Univer- 
sity, Ithaca, New York (Received October 24, 
1958). 


A calculation is made of the fraction of the 1” 
mesonic decays of the hypertriton ,H* which 
yield the two final products He* and 7~. This 
fraction is a function of the spin of ,H* and of 
the ratio p/s of the amplitudes for decay of the 
free A via the s- and p-wave channels. The re- 
sults are compatible with the present experimen- 
tal data. They indicate that probably p/s< 1, 
and that the spin of ,H* is 1/2, which implies 
that the singlet A-nucleon interaction is more 
attractive than the triplet. These results are in 
agreement with those of a previous calculation 
by Dalitz for ,H*. 
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REPRESENTATIONS OF THE DIRAC EQUATION. 
S. K. Bose, A. Gamba, and E. C. G. Sudarshan’ 
Department of Physics and Astronomy, Univer. 
sity of Rochester, Rochester, New York (Re- 
ceived October 20, 1958). 









A representation of the Dirac equation which 
displays its extreme relativistic properties is 
discussed. The “spin” appears naturally related 
to the “position.” 


“present address: Jefferson Physical Laboratory, 
Harvard University, Cambridge, Massachusetts, 


MESON-MESON SCATTERING TERM AND 
STATIC MODEL IN Ps-Ps MESON THEORY. 
M. Sugawara, Purdue University, Lafayette, 
Indiana, and A. Kanazawa, Hokkaido University, 
Sapporo, Japan (Received March 17, 1958). 


A static Hamiltonian has been constructed 
from the renormalizable covariant pion field 
theory in a manner similar to that recently pro- 
posed by one of the present authors. Improve- 
ments are that the meson-meson scattering term 
is included besides the ps-ps coupling term, 
while the latter alone was assumed previously, 
and that an argument is presented to show that 
the Foldy transformation is the unique one ge- 
nerating a valid static Hamiltonian, though it was 
left undetermined before. The resulting static 
Hamiltonian has been analyzed, for the cases of 
low-energy S- and P-wave pion-nucleon scatter- 
ing and threshold gamma-pion production, in 
terms of the one-meson approximation of the 
Chew -Low-Wick formalism without recourse to 
perturbation expansion. It is found in particular 
that the meson-meson scattering term modifies 
the Chew- Low effective range plot of the 6,,- 
phase shift, making the renormalized P-wave 
coupling constant smaller than the customary 
plot gives, for a positive coefficient of the 
meson-meson scattering term in the Hamiltonian 
Empirical values of the coupling constant deter- 
mined through the Chew- Low plot and threshold 
gamma-pion production are shown to be inter- 
pretable in terms of the renormalized P-wave 
coupling constant of ~ 0.08 and the meson-mes0 
scattering term with coefficient of = +7 (h=c=!). 
The present treatment of threshold gamma-pio 
production is criticized in comparison with the 
relativistic dispersion relation. General feature 
of the static model resulting from the ps-ps mé& 
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m theory are summarized in the final section, 
together with the conclusions obtained. The ef- 
fects of strange particles and of renormalization 
have been neglected. 


SOLUTION OF MAXWELL’S EQUATIONS IN 
TERMS OF A SPINOR NOTATION: DIRECT AND 
INVERSE PROBLEM. H. E. Moses, Institute of 
Mathematical Sciences, New York University, 
New York, New York (Received September 15, 
1958). 


Maxwell’s equations for fields with sources in 
media in which the dielectric constant and per - 
meability are unity are written in terms of a 
spinor notation which resembles the one used for 
Dirac’s equation for the electron. One can intro- 





duce Green’s functions and expansions in terms 
of complete sets of orthogonal functions, analo- 
gous to those used in the quantum theory of the 
electron, to solve Maxwell’s equations in more 
compact form than in terms of the convéntional 
vector notation. 
In addition, the new notation enables us to solve 
in a simple way an “inverse radiation problem” 
which we describe as follows: Consider at time 
t<0 the electromagnetic field to be zero. At time 
t=0 sources are turned on and then later turned 
off. The electromagnetic field, which results 
after this process has been completed, will be a 
radiation field. We can solve the problem of find- 
ing the nature of the sources which will lead to 
a prescribed final radiation field. It is shown 
that, in general, the sources are not unique but 
additional conditions can be given which will make 
them so. 









































